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ABSTRACT 


A general mesh generation code (MESHGEN) and finite ele- 
ment flow solver (TURBO) for calculating the flow development 
through axial turbomachines are fully documented. The finite 
element approach followed Hirsch and Warzee. Excellent re- 
sults were obtained for the NASA Task-1 compressor operating 
with subsonic flow conditions. Construction of the code will 
allow straightforward extension to transonic flows, turbine 
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ie Na ODUCT LON 


Pee OLCALEMENT OF TASK 

The original task of this research project was to con- 
tinue the work of Macchi [Ref. 1] and Cirone [Ref. 2] in the 
development of a turbine prediction computer code for the 
Turbopropulsion Laboratory at the Naval Postgraduate School. 
An analysis of the referenced code by Ferguson [Ref. 3] in- 
dicated that a significant amount of work remained to be 
done in order to make the program operational. In the au- 
thor's Opinion the task could be better accomplished through 
the use of a different solution technique. After additional 
study and review of similar work [Refs. 4, 5, 6 and 7] it 
was decided that a finite element approach to the problem 
would be adopted. A program developed by Gavito [Ref. 3], 
which followed the work of Hirsch and Warzee [Ref. 4], was 
selected as the basis for development of the computer code 
Gescribed in the sections that follow. However, Gavito's 
program was formulated as a compressor performance precdic- 
tion which, as it was reported, had not given results simi- 
lar to those obtained by Hirsch and Warzee. Thus the first 
goal of the project became the development of an axial com- 
pressor prediction code that could produce results compara- 
ble to those published by Hirsch and Warzee. The second 


goal was to revise and document the program so that its 


i 





application to any compressor and its extension to turbine 


analysis could be carried out without excessive difficulty. 


Eee OE SCRIPTION OF THE PROBLEM 

One purpose of conducting a through-flow analysis is to 
predict the performance of a turbomachine under design and 
off-design operating conditions. Through the combination of 
a mathematical model and empirically determined correlations, 
it is possible to provide the engineer with a tool that will 
determine the effects of variations in design parameters 
and analyze the performance of an existing machine. 

The first problem addressed in the formulation of a per- 
formance prediction code is that of expressing the analysis 
in a form that can be accurately and efficiently solved by 
computer methods. Most methods for through-flow calcula- 
tions are basec on the classic work of Wu [Ref. 9] which 
stated that the equations of fluid flow in turbomachines can 
be solved on the intersecting sets of stream surfaces known 
eemeene Sl family and S2 family of stream surfaces (Fig. 1). 
In general the intersection of a Sl surface and a S2 sur- 
face is a twisting line with three dimensional variations. 
However, if an axisymmetric assumption is made, the S2 sur- 
face will lie on a meridional plane and the directional 
derivatives on the $2 surface become the 0( )/dr and 9d( )/93z 
mumeylLindrical coordinates. As shown by Smith [Ref. 10], 


circumferentially averaged equations with an axisymmetric 


16 





flow assumption can be used to a good first approximation 
for the through-flow analysis. 

Three general methods of solving the so-called radial 
equilibrium equation of flows in turbomachines which results 
from the axisymmetric assumption can be found in the litera- 
ture. The first is called the streamline curvature method 
(Refs. 1, 2 and 11]. The method derived its name from the 
fact that the radius of curvature of the streamline is an 
integral part of the formulation. The second, a matrix 
method, applies a finite differences technique to the radial 
equilibrium equation, normally after it has been reduced 
to a Poisson form [Refs. i2 and 13]. The third method is 
the finite element method which was used in the present 
work. 

In the mid-1970's, Hirsch and Warzee [Ref. 4] first re- 
ported the application of the finite element method to solu- 
tion of the axisymmetric radial equilibrium eguation. They 
applied the finite element technique to solve the equation 
expressed in quasi-~harmonic form in terms of the stress 
function. They published extensive comparisons of the pre- 
dictions obtained using their method with measurements ob- 
tained on several machines under various operating conditions. 
In general, the method produced excellent results for com- 
pressors and turbines of single and multi-stace configura- 


tions. It was this cemonstrated ability of the method over 
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such a wide range of parameters that let to its selection 
for use in the present project. 

In the sections which follow, the development of vro- 
grams MESHGEN and TURBO, which are based on the work of 
Hirsch and Warzee [Ref. 4], is documented. Comparisons are 
given of the results obtained when the program was applied 
to analyze the NASA Task-1 compressor with results obtained 


by the referenced authors. 
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Il. MATHEMATICAL MODEL 


The equation of motion for a fluid has the general form 
(Vavra [Ref. 14]) 


(9v/ot) + (V-¥)V = -Vp/o + a ~ (gz) (1) 


Using the vector identity 


Conk = vig memy x Vx V) (2) 
Eq. (1) can be written as 
aV/at + V(¥?/2 12) (=e V D7 e+ oe + (VxVxV) (3) 


The first law of thermodynamics for a fluid particle can be 


written as 
Tds = dh - dp/o (4) 


which, along an elemental path length dr in a fluid field 


implies that 
T(dr-Vs) = (dr-V)h - (dr-V)p/o (4a) 
Og 
dr - (Vh-TVs-Vp/o) = 0 (5) 


Since dr is not equal to zero in general, in a homogeneous 


Fluid flow the first law may be expressed as 
Vn = TVs = Vp/o = 0 (6) 


Biestltuting Eq. (6) into Eq. (3) yields 


ee 





s7/dt + V(h+V*/2+ gz) = TVs + £. + UxVx¥ (7) 
For steady, inviscid flow Eq. (7) reduces to 
VH = TVs + (VxVxV) (8) 


As shown by Hirsch and Warzee [Ref. 15], Eq. (8) can 
be revised to describe the flow through blade rows by intro- 
ducing a circumferential averaging process and assuming that 
the flow is axisymmetric at the averaged condition. The 


averaged equation can be expressed as 


~(VxVxvV) = TVs - VH + 3 + a (9) 


where PY is the body force representing the action of the 
blades on the flow and Ps represents the dissipative force 
whose work generates the irreversible entropies. The Pa 
forces are considered to be uniformly distributed in the 
tangential direction and proportional to the loss coeffi- 


cients. Equation (9) leads to the following three equations 


i cylindrical coordinates (with 3( )/d@ = Q) 


Ree) (a teN far) > V_ (lav / oz) - (aV_/ar)) = 


Sate eyo cise rk) = DS ae = Date (10a) 

(ye/x) (ote) /dz) + (l/r) (a(rv)/ar) = F, (10b) 
V_((@V/32) - (8V,/3xr)) - (V/r) (3 (rv) /3z) = 

SH/oz2 - T(ds/sz) - F, Gres) 
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mederon (10a) expresses the radial equilibrium of the 
meridional through-flow and it is the governing equation to 
be solved for the velocity components. Equations (10b) and 
(10c) determine the tangential and axial components of the 
forces once Eq. (10a) has been solved. 

For the solution of Eq. (10a) to have physical meaning, 
care must be taken to ensure that continuity 1s satisfied 
throughout the field. In general, the continuity equation 


can be expressed as 
3o/3t + Ve (pv) = 0 (11) 


which for steady, circumferentially averaged flow can be 


written as 
(9/3r) (prbv_} ce (3/32) (prbv,) = 0 (i2) 


where b is a blockage factor describing the reduction in the 
Flow area in the tangential direction due to the presence of 
rotor and stator blades. The tangential blockage is ap- 


proximated by 
b=i1- t/s (12a) 


where t is the blade thickness and s is the blade spacing. 
As will be discussed in the description of subroutine INPUT, 
this factor will be modified to account for the end-wall 
boundary layer contractions. A stream function, w, can be 
introduced and defined so that Eq. (12) 1s automatically 


Satisfied as follows: 


on 





ow/oar orbV, (3a) 


aw/daz ~orbV,. (13b) 


Peeeer the substitution of Eqs. (13a) and (13b) into Eq. 
(10a) one is assured of the implicit satisfaction of con- 
memeaty as the radial equilibrium equation is solved ex- 


peitliy. Equation (10a) may now be written as 


(9/dxr) ((1/orb) (ap/ar)) + (8/32) ((1/prb) (a¥/az)) = 
(1/v,) ((38/az) - T(ds/dr) + 


(Vv /r) (8(rV,) /8r))- FL > Fay (14) 


The equation is written in a slightly different form for 
solution in rotor regions where rothalpy remains constant 


along a streamline. The definition of rothalpy, H,, given by 


H = ho = CEs (15) 


is substituted into Eq. (14) and the terms in brackets on 


the right-hand side become 


[aH,/ ax - T(as/ar) - (Wy/x) (S(rVy)/3r)- 7, L- Fy. | (16) 
Geessignificance of the Fy = and Fy ” terms can be analyzed 
in the following manner. The body force Fi acts in a direc- 


tion normal to the mean blade surface, which for radial 

blading is in the circumferential direction. The term ae 
f 

accounts for the body force component in the radial direc- 


tion that results when blade lean is present. For most 


blading this is a small term that may be neglected. 


Ox pa 





cemilarly, Pay is the contribution of the dissipative 
forces in the radial direction for non-cylindrical stream 
surfaces. This contribution can normally be neglected for 
axial-flow machines, which is the case treated here. (Note 
that the two body force terms are included in the analysis 
for machines in which the magnitudes of these forces are 
Significant.) With these simplifications the radial equi- 


librium equation may be written in the form 


(3/ar) ((1/prb) (ap/ar) ) + (9/32) ((1/prb) (3p/9z)) = 


ey y_) ((oH/dr) - 1(9s/or))+(V, /xr) (9 (rv) /3x) eer) 


with the appropriate modifications for solution in a rotor 


Pegion. 
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ieee eee BoM ENT METHOD 


Pee ENTRODUCTION 

The finite element method is a numerical procedure that 
is particularly well suited to solving problems in continuum 
mechanics, which invariably involve equations expressed in 
differential form. As stated by Cook [Ref. 16], the essence 
of the finite element method is the "piecewise approximation 
of a function >, by means of polynomials, each defined over 
a small region (element) and expressed in terms of nodal 
values of the function." 

In order to understand the finite element method and 
the solution techniques employed in the computer program 
reported herein, one must first have a complete understand- 
maeeor the basic element, the terminology used to describe 
the element, and the relationship between the element and 
the remainder of the solution domain. The complete problem 
is solved in a piecewise manner, in which the sclution of 
the derived governing relationship over the discrete region 
of an element is sought to determine the contribution of the 
element to the overall solution. Figure 2 illustrates the 
Single element as it is used in the present work and the 
nomenclature for the element on what is referred to as the 
"local domain". The number scheme to employ is arbitrary, 


limited only by the requirement that the system remain 
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consistent from element to element. Figure 3 shows a ver- 
tical stack of three elements to show how elements are con- 
nected in what is known as a “global domain". Table 1 lists 
the relationship between the two reference systems, known as 
the connectivity. The connectivity 1S important because the 
solution of a problem over a computational region involves a 
careful summation of the local contributions of each element 
to the global equations. The summation process is tracked 
by the connectivity. Again the global numbering scheme is 
arbitrary, influenced mainly by considerations of computer 
storage and computational efficiency. 

The key concept to be grasped is that the finite element 
method is a series of local solutions that are coupled to- 
gether through the connectivity relationships to form a so- 
lution for the entire computational domain. A more detailed 
description of the finite element method is contained in 


Beco. 16 through 18. 


Peeeeere Pr LICATION OF THE WEIGHTED RESIDUAL PROCESS 

A standard weighted residuals process was used to trans- 
form Eq. (17) into a form that can be solved by numerical 
techniques. As a first step, the equation was written ina 


more compact form as 


(d/or) [k(dv/or)] + (38/dz) [k(ov/dz)] +£ = 0 (18) 
where 


Klee =e erp) (59) 


BS 





mic, z) = (l/V){Tlds/dr) - dH/ar t+ (V/r)(a(rV) fdr) | (20) 
with the boundary conditions 
k(dw/dan) + a, (b- vo) = 0 (Zz) 


on the associated exterior surface S. Equations (18) and 


(21) may be rewritten as 
(Peay toy om) [keow/ or) }] + (0/dz) [k(ov/oz)] + £} = 0 (22) 
in the volume, V, and 


iste wom) 6+ a, (b - v4) ] = 0 (2.3) 


~?- 


on the surface, S. An approximation, w(r,z), of the unknown 
solution is searched for such that the corresponding weighted 
meendual, R, is equal to zero. Analytically this is ex- 


pressed as 


R = {[W(r,z) blyalar 4) mach a Has piza le eRe (tees) acl Se (24) 
v S : 


where W(r,z) is the (Known) weight function and Ry and Ro 


are the so-called "residuals”™ in the volume and on the sur- 


face, respectively. As the sum of Ry, and R, approaches 


o) 
zero, the approximation, b, approaches the exact solution, 
WwW, with Ry and Ro defined to be identically zero if wy = . 
By defining R,, to be equal to the left-hand side of Eq. 

(22) and R, to be equal to tEhewlett=hane side of Eq. (23), 


pq. (24) can be written as 


f-W(r,z) ((9/3r) {k(3w/dr) } + (38/32) {k(0W/5z) } + EJ] 2mdg 
Q 
+ {Wk (3W/3n)27dC = 0 (25) 
C 
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iiieegration of Ene £lrst term of Eq. (25) by parts yields 


f(k{ (op/or) (aW/ar) + (3~/3z) (oW/dz)} -WE] dQ = 0 (26) 


J 


if w is selected to equal Vo along the corresponding part of 
the boundary. The second term of Eq. (26) reduces to zero 
through the proper application of the boundary conditions. 
The boundary conditions must be satisfied in different 
ways for different portions of the boundary. Along the in- 
let where a, = 0 the conditions may be satisfied by speci- 
fying (sW/dn) to be zero or by specifying the nodal values 
Of W if the inlet conditions are conducive to calculating yp 
for each node. Along the shroud and along the hub the value 
of w must be specified as wy = (m/27) at the shroud and wy = 0 
@eeeene hub. For nodes at the exit plane, the condition that 


(dW/3n) = 0 is required. 


See ee PLLCATION OF THE FINITE fSLEMENT METHOD 

The first step is to discretize the region into sub- 
regions or elements. Within each element the unknown stream 
Function, Ww, and the coordinates r and z are assumed to have 


the following polynomial variations: 


n 
wi{r,z) = d VN, (8,0) (27a) 
n 
r= ) rN, (2,n) (27b) 
a0 
n 
Zz = } ZN, (&,n) (27c) 
Be 


Pad; 





where n = number of nodes in the element 


Ny = the shape or (trial) function for node i 
Ys = value of UW at node 1 
a= value of r at node il 
cea value of z at node i 


Equations (27b) and (27c) imply a geometrical as well as 
functional transformation or mapping, as shown for the pres- 
ent case in Fig. 4. 

The second step in the process is the selection of the 
weight and shape functions. The shape functions are defined 
when the particular type of finite element 1s selected for 
the computational grid [Ref. 16]. The eight-noded quadri- 
lateral was used in the present program and its associated 
shape functions were entered in a subroutine. The weight 
function is independent of the shape function and may be 
chosen at the discretion of the individual. In the present 
case the standard Galerkin technique was employed and there- 
fore, the weight functions were defined as being equal to 


the shape functions, so that 
Wa 2) = Ne 7)) (28) 
Equation (26) may now be expressed in the following form: 
n n n 
fis[iosy 321} cox, /01+ 009,732) | 6,752] 6m, 52)] 


iz 
E 


= NE (e,2)} 42 = 0 (29) 
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where { represents the integral over an element. In matrix 
EB 
Moecation this becomes 


[KI {S}° = {£}° (30) 
where 
Kis = emia) fan foe) ON, / ar) + (GN /92) (3N;/92) 140 (31a) 
eS = IN; f(r,z) ao (31b) 
and 
a (aie) 


The summation of the elemental contributions over the entire 
region yields the global system of equations needed to solve 
the problem. In matrix notation the global system of equa- 


tions 1S expressed as 


Perc = ver (32) 
where 
a e 
Use) iS) 2 (33a) 
1 
AY e 
co d (Sh. (33b) 
a e 
eae ae eae (336) 
Ef 7! 
and m = number of elements in the mesh 
— = p. 


ak ab 


zd 





(K] = system's stiffness matrix 

{f} = system's right-hand side vector 
Since K and f. depend on the unknown solution », Eq. (32) 
is a nonlinear differential equation to be solved by an 
iterative procedure. The details of the procedure are con- 


“ened in section V. 


D. NUMERICAL INTEGRATION TECHNIQUE 

In general, problems are analyzed using a coordinate 
system in which the boundary conditions can be written and 
satisfied in the simplest possible way. For problems with 
irregularly shaped boundaries and/or mixed conditions along 
different portions of the boundary, obtaining numerical so- 
lutions in the original coordinate system can be a formida- 
ble task. Very often a scheme must be found to transform 
the derived equations into a new coordinate system that con- 
forms to the requirements of standard numerical techniques. 
Traditional transformaticn techniques tend to be complicated 
exercises in algebra that require extensive reformulation 
for each geometry or type of boundary condition. The power 
emcee CLinite element method is the automatic inclusion of 
a transformation of the geometry and the function tc a rec- 
tangular domain where a variety of integration techniques 
may be employed. This can be seen in Fig. 4, which illus- 
trates what is implied by Eqs. (27a), (27b), and (27c). The 


method can handle extremely complicated boundary conditions 
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and shapes with ease and is limited only by the type of ele- 
ment selected by the individual. 

The quadratic properties of the eight-node element allows 
curved boundaries in the physical domain so lcng as the sec- 
tion of the boundary included within a single element does 
not have a point of inflection. The use of a quadratic ele- 
ment also ensures continuity of the approximated function 
along the elemental boundaries regardless of the direction 
of approach from within the mesh. The specific numerical 
technique used in the program is discussed in the following 
Peocc1on. 

1. Stiffness Matrix Evaluation 

In section C the following relationship was derived 


for the individual elements of the stiffness matrix, [K]: 


ko. = fk(r,z) ((8N./3r) (SN. /3r) + (8N./9z) (3N,/3z) ]dQ2 (31a) 
9 4 5 i 5 oh 


In order to take advantage of well established numerical in- 
tegration techniques, Eq. (3la) must be transformed from the 
(z,r) domain and its irregular elemental boundaries to the 
rectangular (€,n) domain. Equations (27a) through (27c) de- 
wermoe the variation of the function and the (r,z) coordinate 
femmes in the (€,n) plane. In order to transform Eq. (3la) 
it is necessary to determine the relationship of the varia- 
tions of the derivatives in the two domains. These relation- 
ships can be derived in a straightforward manner through the 


use of the chain rule as follows: 


3k 





(ON, /35) (ON, /9Z) (92/96) + (ON. /or) (dr/9d6) (34) 


and 


(ON, /9n) (ON, /d2) (dz/dn) + (dN, /dr) (dr/én) (25) 


Equations (34) and (35) can be combined in matrix form as 


any faz ar] ( “Ni 
0 & 0€& 0€& 02 
= (36) 
anf [az ar} | 2M; 
dn an on or 


Through the selection of the type element to be used in the 
mesh, N; (&,n) is a known function [Ref. 16], which makes 
possible the computation of the left-hand side vector for 

any point within the element boundaries. Similarly, by tak- 
mogmine Appropriate derivatives of Eqs. (2/7b) and (27c) the 
2x2 matrix, known as the Jacobian matrix yale can be deter- 
mined. It follows that the r and z derivatives of the shape/ 
weight functions can be determined for any point of an 


element from 








al 
COL) Govier We etl ig (pce 
yz re We ae 
= (37) 
BL laa ee a 
or an on an 


Pomexamination of Eqs. (27a) through (27c) shows that once 
the derivatives of the shape/weight functions are known for - 


a point then it is a simple procedure to determine the 


ye 








derivatives of any other elemental property that has a value 
specified at the nodes. 

The final relationship that is needed for the trans- 
formation is the relationship between the differential change 
in the coordinate directions of the (r,z) plane and the (é,n) 
plane. As shown by Kaplan [Ref. 19], the required relation- 
saop is 

Slagle 3" ||| sp [elerelig (38) 


Equation (3la) can now be transformed for integration in the 


(€,n) plane to the form 


aed 
ae y Nk ({B]"{B]) |g |dea (39) 
ij K kk i 
-l -1 
where 
an 
1} % 
fe] - (3 Ho 
oN 
an 
and oe ey A rena of) (41) 
The Gauss-Legendre method of numerical integration 
was used to obtain a solution to Eq. (39). It was selected 


because its determined accuracy was easily determined and 
the simple summing procedure used in the solution could be 
efficiently coded. A one dimensional example is used here 


pomeLiustrate the use of the method. 
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The definite integral 


1 
r= f o(édé (42) 
-1 
may be written in the form 
Ie Wio4 + W565 aa W3 > = a Wren (43) 
where $9, = o(§;) 
W, = Gaussian weight fuction for o5 


The values of the points, called abscissas, and their cor- 
responding weighting function values are catalogued for use. 
The number of points to be used is determined by the order 

of the function to be approximated. In general, a polynomial 
of (2n - 1) is integrated exactly by the use of n Gauss 


points. In two dimensions the Gauss method can be written as 


i eel 
t= | f 6(&,n)dédn (44) 
-l -l 


which can be written as 


tH 
beer 


m 
~ 3} WiW5 (Es ,n5) (45) 


Pomadtton (39) can now be written ina form that can be coded 


Senesolution by the computer as 


k.. = ) Jw g 2 NgKq Veena (46) 
K 
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The scheme used in the program 1s a two dimensional, three- 
point Gauss integration. A detailed description of the 
evaluation of [B] (Bq and the method used to obtain a 1s 
contained in the description of subroutine STIFF. 
2. Right-hand Side Vector Evaluation 
In Section C the following relationship for f was 
derived: 


2 Wn PS Ge es arehe (31b) 
E aE 


By uSing the techniques of Section Dl, Eq. (31b) can be re- 


placed by 
e 
cee Won N 2 NE [Ay | (47) 
a8 5k i Fe tg Naney's 


A more detailed description of the specific methods used to 
Pomve Eq. (47) is contained in the description of subroutine 


BEAL. 
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IV. DESCRIPTION OF PROGRAM MESHGEN 


A. MAIN PROGRAM DESCRIPTION 

MESHGEN was developed in order that the required inputs 
for the program TURBO could be generated in a fast, accurate, 
and conceptually correct manner. The program generates an 
eight-node l1soparametric element mesh, computes the related 
connectivity matrix, defines the type of region enclosed by 
each element, computes the tangential blockage factor and an 
estimate of the stream function for each node in the mesh, 
and computes the thermodynamic conditions and velocity at 
the inlet. The inputs required to operate the program are 
the mass flow rate, total temperature, and total pressure at 
the inlet, the blading and machine geometries, RPM, Cp, y, R, 
and scaling constants for the plot of the mesh. The blading 
geometries must be coded in the program as a subroutine that 
uses approximations or design information to express the 
blade variables as functions of radius. The user provides 
the other information in response to interactive prompts. 
The program has two modes of operation, one which generates 
a complete mesh and ail of the associated information and 
another which uses a previously generated mesh to compute the 
changes in specific arrays that result from a change in the 


Malet conditions. 
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the program is Completely general and may be used for 
either single-stage axial compressor or turbine, and, with 
very minor modifications, can be expanded for use with multi- 
Stage machines. The mesh size that can be generated by the 
program is limited only by computer storage considerations. 
To use the program for another machine, the user is required 
to replace subroutine TASK1 with a subroutine that can com- 
pute the tangential blockage factor, b, for the desired 
blading. The functioning of the program and its subroutines 
for both modes of operation is described in the section B. 

The program's algorithm in outline form is as follows: 

Algorithm: 

Determine the value of the appropriate operating condi- 

tions and whether a new mesh is desired (Subroutine 

NETL). 

Obtain the coordinates of the super element corners and 

a description of the division of the super elements into 

the final mesh. Compute the storage allocation parame- 


ters (Subroutine INPUT). 


Compute the (r,z) coordinates for all nodes in the mesh 
(Subroutine TMESH). 


Compute the connectivity relationships for the mesh and 
determine the beginning and ending node numbers for the 
m@eOr andy stator (Subroutine CONEC). 


Compute the array of node numbers where the value of wb 
1s to be specified. Compute an initial estimate of the 
nodal stream function distribution and call subroutine 
FLOFCT to determine the inlet conditions (Subroutine 
myET 2). 


Compute the nodal tangential blockage factor, b, for the 
rotor and stator nodes (Subroutine TASK1). 


Place the computed values in disk storage (Subroutine 
FILGEN). 


BYP 





Plot the generated mesh on the Tektronix 618 terminal for 
mispection (Subroutine MPLOT) . 


END 


Bee SUBROUTINE DESCRIPTIONS 
me 6 SUbEeuUcine INIT) 

Subroutine INIT1 obtains the value of thermodynamic 
variables and plotting parameters that are required for the 
program in either mode of operation. The user is required 
to provide the values of the mass flow (lbm/sec), total tem- 
perature (°R), and total pressure (psia), ratio of specific 
heats (vy), gas constant (ft-lbf/lbm-°R), and the specific 
heat at constant pressure (BTU/lbm-°R). The scaling con- 
stants are convenient values of r and z used to frame the 
plot of the mesh. If a new mesh is not desired, the program 
exits the subroutine. 

Subroutine INIT1 determines if a new mesh is to be 
generated by the response of the user to an interactive gues- 
tion. If a mesh is to be generated, the subroutine obtains 
some preliminary information about the flow regicn. Figure 5 
shows how the user must first divide the flow region into a 
coarse mesh known as super elements, recording the (z,r) co- 
Ordinates of the corner points. The minimum number of super 
elements for a single-stage compressor is five so that the 
three duct regions, the rotor, and the stator may be repre- 
sented. The maximum number cf super elements and the subse- 


quent division into mesh elements is limited only by the 
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storage limitations of the machine. In practice, the maxi- 
mum number of elements is limited by the number of equations 
that may be solved by the program TURBO. It is also limited 
by the fact that only one super element may be used to de- 
scribe a rotor or duct region and that the super elements for 
these regions may only be further subdivided into a single 
column of elements. The latter restrictions are for compati- 
bility with the computation procedures used in the program 
TURBO. A decision must then be reached on what subdivision 
of the super elements will provide a mesh that is sufficient- 
ly fine to yield accurate results efficiently. Once the flow 
region has been divided into super elements and a determina- 
tion as to the total number of rows and columns of mesh 
elements has been made, the user can input the appropriate 
values in response to the prompts provided by the program. 
2. Subroutine INPUT 

Subroutine INPUT uses interactive prompts to obtain 
a description of the turbomachine's flow passage geometry and 
the desired mesh characteristics. The user must provide the 
program with the coordinate values of the super element nodes 
as shown in Fig. 5. The values are entered as nodal pairs on 
a station-by-station basis. The first (z,r) coordinate 
entered is the node on the shroud and the second lies on the 
hub. The program then asks the user to identify the type of 
region enciosed in each super element and into how many 


columns each super element is to be divided. 


aS, 





Enough information is now available for the program 
to compute and store the information required for the program 
TURBO. Subroutine INPUT stores the responses to the inter- 
active prompts, determines the values of the storage location 
pointers, and determines if any stcrage limitation has been 
exceeded. If any storage limitation is exceeded, the sub- 
routine calls the appropriate error subroutine to halt execu- 
tion. The interactive portion of the subroutine is omitted 
if a new mesh is not desired. However, the values of the 
pointers are calculated and storage requirements are evalu- 
ated as before. A listing of the pointers and the corres- 
ponding variables is given in Appendix A. 

ee Subroutine TMESH 

Subroutine TMESH computes the nodal coordinate values 
from the inZormation obtained in subroutines INIT1 and INPUT. 
The subroutine uses linear interpolation in the axial direc- 
tion and a quadratic scheme in the radial direction. The 
radial interpolation scheme maintains the difference in the 
squares of the radius of the nodes as a constant. This al- 
lows the assumption of equal mass fiow between the nodes for 
uniform axial velocity which is used to determine the initial 
estimate of the nodal stream function distribution. The 
nodes are numbered with the assumption that the fiuid flow 
is from left-to-right in the mesh and that the number of 
columns of elements is greater than or equal to the number 


of rows of elements. The node at the inlet shroud is labeled 
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1 and the node at the outlet hub is labeled n for ann node 
mesh. The numbering proceeds on a column-by-column basis 
from top-to-bottom. The total number of rows, columns, and 
mesh elements is displayed to the user. The mesh computa- 
tions are omitted if a new mesh is not created and elemental 
count information is displayed as before. 
ie Solibroutine CONEC 

Subroutine CONEC generates the connectivity matrix 
for the mesh. The connectivity matrix is used to keep track 
of which nodes are assigned to which elements and the ar- 
rangement of the assigned nodes within the element. The con- 
nectivity relationships for a three element stack is shown 
in Fig. 3. Additionally, the subroutine determines the first 
and last nodes of the rotor and stator. 

pee ocubroutine INIT2 

Subroutine INIT2 stores the node numbers for nodes 
where the value of » is specified as a known quantity. The 
array 1s used in the program TURBO to apply the boundary 
conditions. The array is not computed if a new mesh is not 
desired. For either mode of operation, subroutine INIT2 
computes the values of the inlet thermodynamic variables, 
the inlet axial velocity, and an initial estimate of the 
nodal stream function distribution. Subroutine FLOFCT is 
used to calculate the iniet conditions and is described in 
the next section. The initial stream function is computed 


from the boundary conditions at the shroud and hub. Along 
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the shroud, w is specified to be equal to (m/2m) and along 
the hub to be zero. The value of w along the inlet is de- 
termined by a linear interpolation because of the quadratic 
node spacing in the radial direction. The remaining nodal 
values of W are obtained by an assumption that ww 1s a con- 
stant along the streamwise boundaries of the elements. The 
assumption is obviously in error, but it observes the boun- 
dary conditions and provides a reasonable first estimate to 
begin the iteration scheme used in the program TURBO. 
ee subroutine FLOFCT 

Subroutine FLOFCT computes the inlet conditions for 
a passage with a specified geometry, mass flow rate, total 
temperature, total pressure, and an assumed uniform inlet 
velocity. The method followed is the "total flow function" 
formulation proposed by Shreeve [{Ref. 20]. The total flow 
function is defined as the ratio of the mass flux to the 
Pmmeing Or stagnation mass flux. The following definitions 


and equations are required for the method: 








ee 0.5 
Vi. = ({2H] (48) 
oe V/V, (49) 
iy ee Xe (50) 
% 
5 (y-1) 
p/P, ee (le (S10) 
1 
(y-1) 
p/ po, = (1 - X”) (a2) 





Mmeem the definition of the total flow function, 6, it follows 
that 
y-l 


6 = pV/o,V, = X(1 - X*) (53) 


The value of » at the inlet can be calculated at the inlet 


from the assumed uniform conditions by the expression 
>, = m/ (0, V,A) (54) 


The value of X at the inlet is found through the following 


Newtonian iteration: 


4 = m/ (0, VA) (54) 
Pesume X = 0.1 to assure the selection of the subsonic root. 
aes 
2 71 
Calculate: @= X(1 - xX) (33) 
one ay eve oe = yl = x7) 1) 
erst |o, SK Te 


Ti the test fails then calculate 


= eee (>, - >) (do/ax) 


and recalculate $? until convergence is reached. Once con- 
vergence is reached the inlet conditions are computed by 
eeaeions (50) through (52). 
fe SUubrOurmnesTASK 1] 
Subroutine TASK1 computes the nodal tangential block- 
age factor for the blading of the NASA TASK1 transonic com- 


pressor. The value of the blockage factor is determined by 
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b=1-t/s (12a) 


The values of t and s are obtained by approximations to the 
known blading geometry that are expressed as functions of 
radius. The maximum thickness of the blade 1s artificially 
defined to be at the mid-point of the chordline to ensure 
that the factor 1s accounted for in the calculations in the 
program TURBO. This artificiality could easily be removed 
through a modification to the axial interpolation scheme used 
in rotor and stator super elements. Subroutine TASK1 is the 
only machine dependent subroutine in use in the program and 
would need to be replaced with an appropriate substitute in 
order for the program to be used on another machine. The use 
of the subroutine is omitted if the user does not desire a 
new mesh. 
8. Subroutine FILGEN 

Subroutine FILGEN places the computed mesh parameters 
on disk storage for use in the program TURBO. If the limited 
mode of operation was selected by the user, the subroutine 
Suey Updates the values of the parameters that change for a 
new inlet condition. A listing of the output parameters and 
their corresponding storage location 1S given in Appendix B. 

9. Subroutine MPLOT 

Subroutine MPLOT provides an on-line plot of the com- 
puted mesh on the Tektonix 618 graphics terminal. Figure 6 
shows the 63 element, 222 node mesh used in the computations 


of the test cases. The subroutine displayed the mesh through 
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m@merect calls to the subroutines of the library plotting 
package, GRAFF. The subroutine's algorithm is as follows: 
Algorithm: 
Form two Real*4 arrays from the information in the r co- 
ordinate and the z coordinate arrays for plotting 


Gompacibility. 


Sort the arrays and plot the streamwise boundaries of the 
mesh elements. 


Sort the arrays and plot the transverse boundaries of the 
mesh elements. 


END 
me ~Subroutine ERR1I 
Subroutine ERR1 is called by subroutine INPUT if the 
storage limitation for Real*8 variables has been exceeded. 
The subreutine displays the amount by which the limitation 
was exceeded and terminates the program's execution. The 
user response would be to increase the value of LIMR if pos- 
Sible or reduce the size of the mesh. 
mie Subroutine ERR2 
Subroutine ERR2 is called by subroutine INPUT if the 
storage limitation for Real*4 variables has been exceeded. 
The subroutine displays the amount by which the limitation 
was exceeded and terminates the program's execution. The 
user response would be to increase the value of LIM4 if pos- 
Sible or reduce the size of the mesh. 
M2. Subroutine ERR3 
Subroutine ERR3 is called by subroutine INPUT if the 


storage limitation for Integer*4 variables has been exceeded. 
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The subroutine displays the amount by which the Limitation 
was exceeded and terminates the program's execution. The 
user response would be to increase the value of LIMI if pos- 


Sible or reduce the size of the mesh, 
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VV. DESSRIPTION OF (RBROGRAM TURBO 


Program TURBO solves the quasi-harmonic stream function 
radial equilibrium equation for flow in an axial compressor. 
The program uses the information computed by the program 
MESHGEN to calculate the desired thermodynamic information 
at all nodal points and displays selected values on a graph- 


wes terminal for inspection. 


A. MAIN PROGRAM DESCRIPTION 


The program obtains a solution of the equation 
poi (£) (32) 


An iterative scheme was adopted for this nonlinear problem, 
whereby an estimate of the stream function distribution is 
used to calculate values of the velocity components and 
thermodynamic variables at the nodes of the mesh. These 
computed values are then substituted into Eq. (32) and a new 
value of the stream function distribution is calculated. The 
estimate of the distribution is compared to the calculated 
value to determine if the solution has reached convergence 


Geeerding to the following criterion: 


n eects 
rr = tah 
an = ilar IL (47) 
—_ 
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estimate of W at node i 


= 
oi 
D 
K 
D 
< 
| 


~ n+] 
oe 


solution for » at node i 


T£ the maximum difference for all nodes 1s less than ¢€, the 
procedure is terminated. If the maximum difference exceeds 
some specified value of ¢, the new estimate of to be used 


for the next iteration is determined using 


a 


y ntl - v2 - af50t 7 vi] (48) 


where ve = new estimate of ) for the next iteration 
a = under relaxation factor required for conver- 
gence because of the strong nonlinear proper- 
tes Gt Eq. (52) . 
The process of constructing the inputs required for Eq. 
(32) is repeated until convergence is obtained. The details 
of calculating the inputs and constructing the stiffness 
matrix and the right-hand side vector are contained in de- 
scriptions of the program's subroutines. 
mee program's algorithm follows in outline form: 


Obtain the computational constants (SUBROUTINE RDATA). 


Determine the values of the pointers used to partition the 
meerage arrays (SUBROUTINE INIT1). 


Set the initial values for all storage locations to 0.0 or 
QO as appropriate (SUBROUTINE ZEROI). 


Recall from storage the externally computed input values 
mi@einitialize the inlet conditions (SUBROUTINE INPUT). 


Calculate a velocity and thermodynamic variable distribu- 


tion based on the assumed stream function distribution and 
the inlet conditions (SUBROUTINE DIST). 
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From the distributions obtained in DIST, calculate the 
right-hand side vector {f} (SUBROUTINE FCAL). 


Using the density and blockage factor distributions, form 
the stiffness matrix [K] (SUBROUTINE STIFF). 


Solve the system of linear equations to obtain a new stream 
mimetlon distribution (SUBROUTINE DSIMQ). 


Place the solution vector in its proper storage location 
(SUBROUTINE REPLA). 


Compare the original stream function distribution to the 
solution vector to determine the maximum difference in the 
SeestEcibutions for all nodes (SUBROUTINE TEST). 

Determine if the convergence criterion has been satisfied. 
If convergence has not been obtained, perform the relaxa- 
tion iteration to update the estimate of the stream func- 
tion distribution (SUBROUTINE RELAX), prepare for another 
@yele (SUBROUTINE NOCON), and then return to SUBROUTINE 
Por tor further calculations. 

If the convergence criterion has been satisfied, print the 
results (SUBROUTINE OUTPUT) and display selected informa- 
tion on the graphics terminal (SUBROUTINE MPLOT). 


END 


Seo UBROUTINE DESCRIPTIONS 
Sections 1 through 23 provide a detailed description of 
the subroutines of program TURBO. 
1. Subroutine RDATA 
Subroutine RDATA is used to store the following com- 
putational constants: 
(a) Logical Input/Output variable NREAD and NWRITE. 
(5) Relaxation factor. 


(c) Limits for the storage arrays. 


49 





(d}) Three-point Gaussian abscissas and weighting values. 
(e) Constants used for conversions between different units. 
De Mfebbepaenoicmigts UNS RAM 
Subroutine INIT1 determines the values of the point- 
ers used to partition the Real*8, Real*4, and Integer*4 ar- 
rays and determines whether the storage limitations for any 
of the arrays has been exceeded. If the storage limitations 
have been exceeded the subroutine will halt execution by 
calling the appropriate error subroutine. A listing of the 
pointers and their corresponding array names is contained in 
Appendix C. 
Be ‘subroutine ZEROL 
Subroutine ZEROI sets the initial value of all arrays 
equal to 0.0 or 0 as appropriate. 
4. Subroutine INPUT 
Subroutine INPUT retrieves required input information 
meomeits Corresponding disk storage location. The informa- 
tion must be placed in storage before running program TURBO. 
The usual method of generating the information and placing 
1t in storage is through the use of the program MESHGEN. 
Subroutine INPUT also initializes the inlet condi- 
tions to their proper values, and modifies the nodal blockage 
factors to account for end-wall boundary layer effects. No 
attempt was made to include a global method for calculating 
the blockage factors; rather, a method similar to the one 


used by Hirsch and Warzee [Ref. 4] was used. The full method 
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used by Hirsch and Warzee was to artificially reduce the 
size of the flow passage by reducing the boundaries of the 
mesh, followed by the application of a general blockage fac- 
tor to the nodes of the mesh. In the program TURBO, no mesh 
modifications are made. The procedure followed was to apply 
a general blockage factor to all nodes, followed by the ap- 
plication of an additional blockage factor to nodes in the 
outer elements in the rotor, stator, and the passage in be- 
tween. Though reasonable results were obtained by this 
method, the handling of the end-wall boundary lavers remains 
the most obvious weakness in the code. This is addressed 
specifically in section VII. 
Pee Subroutine Dist 

Subroutine Dis! calculates the distributions of 
velocity, density, temperature, pressure, fluid flow angles, 
entropy, and enthalpy using the known blade and machine 
geometry, inlet ccnditions, and the assumed distribution of 
the stream function. Properties of nodes at the mid-line 
of the rotor or stator blades were assumed to have a value 
equal to the average of the inlet and exit conditions of the 
blade. The elemental calculations are accomplished through 
mae control of subroutines SLINE, DUCT, ROTO, and STAT. 

The following is the subroutine DIST's algorithm in 
Seaeline form: 

For each element in the mesh. 


Determine type element for appropriate computations. 


Spl 





Duct Elements 
Ber each Node at Stations 2 and 3; 


Determine the location of the streamline and thermo- 
dynamic conditions at station one (Subroutine SLINE). 


(If the element is along the machine exit plane, en- 
sure that (3/3z) = 0.) 


Compute the thermodynamic conditions (Subroutine DUCT). 


Assign the appropriate values to the proper storage 
Tecation. 


Rotor Elements 

For each node at stations Il: 
Determine the location of the streamline at station 3 
and the 3~/3z and the d~/dr at stations 1 and 3 (Sub- 
BOleine SLENE) . 


Determine the inlet and outlet relative flow angles and 
the outlet absolute flow angle. 


Compute the total-to-total pressure ratio and the adia- 
batic efficiency for the streamline (Subroutine ROTO). 


Assign the appropriate values to the proper storage 
ligeae lon. 


For each node at stations 2: 


Determine the location of the streamline and the ther- 
modynamic conditions and the 3/3z and the d/sr at 
Stations 1 (Subroutine SLINE). 


Determine the location of the streamline and the dy/3z 
and the 3W/dr at station 3 (Subroutine SLINE). 


Determine the inlet and outlet relative flow angles and 
the outlet absolute flow angle and the relative devia- 
tion angle. 


Compute the thermodynamic conditions at station 3 (Sub- 
routine ROTO). 


Compute the value of all properties for the node as 


peling the average of the values at station 1 and sta- 
wahoo Ces ae 
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Assign the appropriate values to the proper storage 
LOecat1On. 


For each node at stations 3: 
Determine the location of the streamline and thermo- 
dynamic conditions at station 1 and the value of 3)/3z 
and ot/dr at Stations l and 3 (Subroutine SLINE). 
Determine the inlet and outlet relative flow angles and 
the outlet absolute flow angle and compute the thermo- 


dynamic conditions (Subroutine ROTO). 


Assign the appropriate values to the proper storage 
IKsyoy=henkes a= 


Stator Elements 

Miiemstatornsalqorithm 1s the same as the rotor algorithm 

except the outlet absolute flow angle is the only angle 

calculated. The inlet absolute flow angle is determined 

through interpolation. 

oe Subroutine SLINE 

In order to understand the functioning of this sub- 

routine and others to follow, one must refer to the nomen- 
clature used to describe the eight-node element. Figures 2 
and 3 show the nomenclature clearly and Table 1 demonstrates 
the connectivity. All of the calculations in the program for 
the distributions of velocity, flow angles, and thermody- 
namic properties are founded on the assumption that the 
points in question lie on the same stream surface. Thus the 
objective of the subroutine is to obtain the location of a 
given value of the stream function at a specified station in 


the flow region. The location of the streamline is required 


in order to compute the variables used in [K] and {f£}. 


aye 





Titouga tnesdpolicatrion Of the boundary conditions, 
the nodes along the shroud are defined to lie on one stream 
surface and the nodes along the hub are on another. It is 
possible for all other nodes in the mesh to be on different 
stream surfaces. For these nodes an interpolation scheme 
must be followed to find the (§,n) coordinates of a specified 
stream surface at a given station in the mesh, 

The solution sequence that the program follows starts 
at the top element of the first column of elements in the 
mesh and solves the thermodynamic and velocity conditions 
for all nodes in the element using the assumed stream func- 
tion distribution and the specified inlet conditions. When 
the calculations for the first element are complete, the 
program continues down the column until the calculations are 
complete for the element along the hub. The program then 
sequences to the top element for the next column and con- 
cinues until the calculations are complete for the last ele- 
ment in the mesh. In this sequence it is always possible 
to calculate the conditions at station 1 of an element for 
any interim nodal values of the stream function. 

The process will be described by way of an example 


for one node as shown in Fig. 3. 


Example: Find the (§,n) coordinates of the streamline that 


passes through node 7 of element 1, Node (1,7). 


@eom the connectivity relationships it is known that 
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(Node (1,7)) = (14) 


The value of (14) is known and the search is begun to find 
two nodal vaules of &— at station 1 that bracket the desired 
value, (14). The program first tests to see of (Node (1,7)) 
is greater than (Node (1,5)). In this case it is not and the 
program would automatically shift and test to see of 

(Node (1,7)) is greater than (Node (2,5)). In this example 


the value is larger and the same test would be applied to 


(node (2,4)). Again the answer would be true. The program 
would then test to see if the value of (Node (2,3)) were 
larger than (Node (1,7)). The answer being true would sig- 


nal the program that the location of the streamline had been 
bracketed and a half-interval technique would be applied to 
memcdeche location. As shown in Fig. 4, the value of &— for 
mmeelocations along station 1 is -l. This fact is important 
for two reasons. One, with € known the program is only re- 
quired to iterate on to obtain convergence. Two, the 
Kronecker delta property of the shape functions means that 
only the shape functions at station 1 have nonzero values 
[Ref. 16]. For the half interval method, the program uses 
the average n of the most recent bracketing as its estimate 
for n. In this example the first estimate of n is equal to 


0.5 and the solution for wt at (-1., 0.5) can be written as 


ma 1.,0.5) = N3(-1.,0.5) 0, uE Nygt(-1.,0.5)0, + No(-l.,0-5) 0. 





The solution is then compared to w (Node (1.7)) to determine 
if the difference is less than some «. If the difference 
exceeds ¢«, the new estimate for n becomes 0.25 or 0.75 de- 
pending on whether the solution is larger than or less than 
the value of w(Node (1,7)). The process is continued until 
convergence is reached. Once € and n for the streamline 
location at station 1 are known, all of the properties for 
station 1 can be determined. (The same method is used to 
find the location of the streamline at station 3 for rotor 
and stator elements.) Having determined the coordinates 
Of streamlines at all desired locations it is possible to 
calculate the required (ON, /ox) and (ON, /32). The computed 
inlet and exit coordinates and conditions are then passed 
to subroutine DIST for use in subroutines DUCT, ROTO, and 
STAT aS appropriate for the calculation of the conditions 
at Node (1,7). 
The following is the subroutine's algorithm in out- 

meme LOrm: 

Duct Element 


If the node being investigated is on station one, exit the 
sebroutine. 


For stations two and three, determine the streamline co- 
ordinates and the thermodynamic conditions at station one 
and compute the awW/3z and dt/dr at the node. Exit the 
subroutine. 


Rotor/Stator Element 


If the node being investigated is on station one, set all 
inlet thermodynamic variables equal to the corresponding 
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nodal value and compute the dy/dz and the dy/dr at station 


mpamee. Exit the subroutine. 


For station two, determine the streamline coordinates and 


the thermodynamic conditions at station one 


and compute the 


Streamline location and the dW/dz and the dwW/dr at station 


three. Exit the subroutine. 

For station three, determine the streamline 

the thermodynamic conditions at station one 

dW/oz and the dv/or at station three. Exit 
me «6Subroutine DUCT 


Subroutine DUCT determines the values 


pressure, and density for the elemental nodes 


coordinates and 
and compute the 
the subroutine. 


of temperature, 


at stations two 


and three. An iterative procedure is used with the knowledge 


that angular momentum is a constant in a duct. 


estimate of the velocity at station 1 is made 


The initial 


uSing the com- 


puted values of ov/dz, dv/dor, r, and b at the node (Subrou- 


tine SLINE) and by choosing the estimate of the density at 


the node to be equal to the density at station 1. The fol- 


lowing sequence of calculations is repeated until convergence 


on a value of the exit velocity: 


O25 
cae ra pues 
Ving 7 (U/lporgb5)) * [(ap/az,)* + (3¥/ar,)*] 
= -1 
A> = tan f(r Vy tan 4) /(rgVi9) | 
ue wae 2 eee : 
T. = Tr (or= 1 ju/ Dex (v.49 {1 + tan ay) Ay RG) 
= Pry (To /Tyy) ** (¥/y-2) 
a P,/(RT,) 
5 ) Oee5 
an = (1/(p5r,b,)) x | (3u/32), ~ (op/or), | 
Test if iV - V < oe 


m2n m2 | 


a7 





Hae total conditions are then calculated from the static 
conditions and the computed velocity. 
Bey SUDrOutIne ROTO 

Subreutine ROTO calculates the change in the relative 
flow angles, the velocity, and the thermodynamic properties 
along a streamline across a (compressor) rotor element. The 
program uses the conditions at station one and the location 
of the streamline and the partial derivatives of the stream 
function at station 3, all of which were calculated in sub- 
routine SLINE. The relationships in subroutine ROTO are 
derived from cascade correlations and known property rela- 
tionships for a streamline in a rotor. 

The first step is the calculation of the inlet and 
Seema celative flow angles. For inlet flow without swirl, 


the relative inlet flow angle can be calculated using 


Ere | 
oT | enyeg9 eee 


The incidence angle is calculated using the known blade 


geometry and inlet angle, since 


or 
where Ky 1s the angle formed between the tangent to the blade 
Ghordline and the axial direction. In order to calculate the 


exit relative flow angle one must determine the deviation 


angle, 4. The program uses the correlations and equations 
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derived by NASA [Ref. 21]. The specific sequence of equa- 


evens, using the notation of Ref. 21, is as follows: 


~ 
1 
te 


= £(M,r) 


ke 
l 
A 


(K + no 


i’su'*o) 10 


(3 (gale eae (a 


Oe eS ou. 0) 110 oe a yp 


Sa 7 Sop = F(M,,4) 


eer sop toe ~ ® 2p) 


Leeg) (0 /di) 5, 


W 
il 
W 
t 


i) -i + 6 


The expressions used to approximate the NASA correlation 
curves were those obtained by Crouse of NASA lewis and were 
Peevided to the author by Okiishi [Ref. 22]. 

When the inlet conditions and relative flow angles 
are known it is possible to determine the conditions at the 
rotor element exit. The initial exit velocity is obtained 
in the same way as in subroutine DUCT. The following se- 
quence of equations is solved iteratively until convergence 


for the exit velocity is reached: 


= ae es | 
a, = tan [ (or, Ving tan 85) /Vn2 | 
D=1- (W,/W,) - (rc, Wo - roW 5) / (2 WX) 
Ww = curve fit to w(D,cos 8,0) 
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cee m 2 2 E: Pa 
Tey = T,, + w (x, ry )/constant 
Poy = Pry (Tey Ty) * Y/y-2D) 
2a, Pry (Tay/ Tp) * y/¥-1) 
Pro = Pry ~ ¥fPpy 7 Py) 
W, = Vin /Cos 8, 
- 2 
T, = Toy - W, /constant 
P, = Poo (T,/Tp,)* (y/y-1) 
P54 = P,/(RT,) 

2 2719: 
ee = (1/ (05155) * | (90/32), > (sp/or), | 
Mest. it Vv _ -V | se 

min m2 


When convergence is reached the total conditions are calcu- 
lated from the static conditions and the computed velocity. 


The value of the entropy change is calculated by: 


Ga= Rk In (P,4/ ) * constant 


Pri 


9. Subroutine STAT 
Subroutine STAT calculates the stator element exit 
absolute flow angle and thermodynamic conditions using the 
knowledge that the total enthalpy is a constant across the 
Stator. The initial estimate of the exit velocity 1s ob- 
tained in the same way as in subroutine DUCT. The following 
sequence of equations is used until convergence for the exit 


velocity is reached: 
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G25 
_ 2 2 | 
We = | Vint (tar a4) 
D=1- (V,/V,) + (r5V4 ~ roV5)/(2o0Vr) 
jem setave fic £O w(Drcos 6,0) 
P = Pp 


rl 7 W(Pry - P,) 


T, = Tp 7 (y-1) /2* (V_4*(1 + tan® 5))/(y RG) 


es Pir (T5/Tpo) ** (¥/y-1) 


O05 = P,/(RT,) 
5 5 U2 
Vmin = (1/legryb,))* [(3¥/92) 57 + (30/34) 4 | 
Test if |v... - V | < € 
min m2 


When convergence is reached the total conditions are calcu- 
lated from the static conditions and the computed velocity. 
The value of the entropy change is calculated by the same 
method used by sukroutine ROTO. 
yO. Subroutine FCAL 

Subroutine FCAL uses the previously computed distri- 
butions of total temperature, entropy, enthalpy, axial velo- 
city, and tangential velocity to compute the right-hand side 
vector for the global system of equations. In the absolute 


frame of reference, f(r,z) can be expressed in the form 
me,Z) = eT) Seale otom Eve E) to (ev ) 7 or) (20) 


The value of f(r,z) within an element is determined using 


the following relationships: 
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m2 ,2) = 5 N,(S,n)T, gl Car gy | j N;,(é,n) 8, 
L a 
n n 

Sa 27) } N. (En) Ss; V(z,r) = ) No (EnV, 
1p n 

AZ) d N.(E,n)V,; a= d N.(€,n)xr; 


where the value of N, (&,n) 1s determined by the value of € 
and n for a specified Gaussian integration point within the 
element. The required partial derivatives are found in a 
Simple and direct way. To illustrate, the (csH/3r) is de- 


rived as follows: 


nh Yi 
(QH/3r) =} (9N,/ar)H, + } N, (3H, /3r) 
a i 


and since H; 1s a constant then 
n 
IH/ar = d (ON, /3r)H, 


where (ON, /9or) Peetound by oq... (3 7); “and H. is the appro- 
priate nodal value of H. The radial variations in entropy 
and angular momentum are found in the same manner. The same 
procedure is followed for the values of rothalpy and relative 
tangential velocity for rotor elements. 

It 1S now possibie to calculate the quantities in 
the expression for the right-hand side vector at a point. 


Ali that remains 1s to apply an integration technique to 
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obtain the value over an element and to assemble the re- 
sulting local contributions into the global equations. As 
Bmewn in section III.D.2, the local contribution for node i 


in an element can be expressed as: 
Ber ) LWW), Mia Ny ep (47) 


In the program Eq. (47) is modified to 


f. = re rN iN ‘a plole 
i in m ¢ i K ok 
where the Gaussian abscissas and the corresponding product 
of the weight functions are grouped into three one-dimensional 
arrays. At the completion of the summing process, the local 
contribution for f has been calculated for each node in the 
element. The global system is then updated by adding the 
local contributions to the global values through the use of 
the connectivity relationships. 
The following is the subroutine's algorithm in out- 
line form: 
Algorithm: 
Iterate for each element in the mesh, 
Iterate for each Gaussian point. 
Find shape functions, |J|, and yeas 


Find Vee i ean a rvs (d9s/dr), and (dH/ér). 


ie 3 


Compute the contributions of the value f at the Gauss 
point to the value of £ at each node of the element. 
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Upon completion of the Gaussian integration, add the 
local contribution to the global system. 


END 
The same algorithm is followed for rotor elements with the 


appropriate substitutions of H, and Woe 


R 
wee rouUDrOuLine STILE 
Subroutine STIFF uses the computed distributions of 
density and blockage factors to form the stiffness matrix 


for the global system of equations. It was shown earlier 


that the contribution to the elemental stiffness is expressed 


as 

ko = inotc , Zia oN ./ dale(oN. / 3%) + (IN./0z) (ON. /9z) J] AQ (31a) 
sling) RE J ae 7 a 

where Katie) =p PED} (ro) 


Again, the elemental properties are considered to have a 


Potynomial variation of the form 


n n n 
e(z,r) = N05, re nae aesand biz .0)) = )N.b, 
nu al a. 
As shown earlier, Eq. (3la) can be converted by the Gauss- 
Legendre method to 
e _ cn ind rE 
Ki > T5W,W Tigh | (1 ares (46) 


The value of Kn LH dererminea from the Gefinitions of o, b, 


and r by the same methods used in subroutine FCAL. The 


e 


evaluation of [{B] [B] is a simple matter to perform. The 
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matrix [B] is simply the column vector ( (ON, /9oz), (ON, /or)}. 
Therefore, ie isl can be written es ee) + 

Bee 22) (ONE Oz)) - The value of (ON, /d2) and (ON, /or) is 
found for all nodes in the element in one step through the 
use of subroutine JACOB. The matrix product can be evaluated 


at a point (z,r) as 


K-15 
Uo 9 


es ae) - (OWN CEA AEN eee), 


By using the same Gaussian weighting scheme used in subrou- 


tine FCAL, Eq. (46) may be written in the form 


© n nn 
e _ 3 
ki = LM, st ae a ae + (3N;/3z) (8N,/32) 1} J| 


The resulting 8x8 elemental matrix is then added to the 
global stiffness matrix through the connectivity relationships. 
Up to this point [K] and {f} have been assembled with- 
out regard to the boundary conditicns except at the exit plane 
where the (d3W/3n) = 0 was enforced explicitly during the pro- 
cedures used by Subroutine DIST. Care must be taken to en- 
Sure that the boundary conditions for the other three segments 
of the boundary are not violated. As shown in section III.B, 
the boundary condition for the nodes along the shroud, along 
the hub and at the inlet plane of the machine is that the 
value of ~ is specified. If the value of ~ is specified at 
these locations the Eq. (32) must be modified so that ») is 


no longer free at these nodes. A standard technique is 
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employed to remove individual equations from a system of 
equations when the degree of freedom represented by the in- 
dividual equations has been removed. 
The following is the subroutine's algorithm in out- 
line form: 
Algorithm: 
Iterate for each element in the mesh. 
Iterate for each Gaussian point. 
Find shape functions, |J|, and iat mee 
Find the value of k at the Gauss point. 


Compute the elemental stiffness matrix. 


Upon completion of the Gaussian integration, add the 
local contribution to the global system. 


Upon completion of the addition of the last element's 
contribution, modify the system of equations to include 
the boundary conditions. 
END 
m2. Subroutine DSIMO 
Subroutine DSIMQ 1s a non-IMSL library, double pre- 
Cision subroutine that solves a set of n simultaneous equa- 


tions of the form 
[A] {x} = {B} 
where [A] 1S an nxn matrix 
{X} and {B} are nxl vectors. 
i. subroutine REPLA 
Subroutine REPLA places the solution vector obtained 


from subroutine DSIMQ into its proper storage location. 


66 





ba. ~SubzOurInNe TEST 
Subroutine TEST determines the maximum difference 
in the assumed nodal distribution of the stream function at 
the beginning of an iteration to the solution of the radial 
equilibrium equation calculated using the assumed. distribu- 
tion. The difference in the distributions at a node is de- 


fined as 
ies he iain’ ~ n+1] 
Diff = (v5 7 Jo, | 


where veo is the assumed value at node i and pe 


1s the 
calculated value for node 1. Convergence is considered to 
be reached when the maximum difference for any node is less 
than a selected reference value, e. 
15. Subroutine RELAX 

Subroutine RELAX performs the relaxation scheme to 
obtain an updated estimate of the stream function distribu- 
tion for the next program iteration. The new estimate for 


stream function distribution is calculated as follows: 
oc) ae n io Cl nn rae 0 
1. = 7+ a, b | (48) 


16. Subroutine NOCON 
Subroutine NOCON prepares the program for the next 
iteration by setting the elements of the right-hand side 


meceor, {f}|, and the stiffness matrix, [K], equal to zero. 
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We ecuoroutine OULPUT 
Subroutine OUTPUT prints the computed nodal values 
of a majority of the velocities and thermodynamic properties. 
A listing of the values that are printed and the correspond- 
ing units is contained in Appendix D. A sample output list- 
ing 1s contained in Appendix G. 
18. Subroutine MPLOT 
Subroutine MPLOT uses the Tektronix 618 terminal to 
make an online graphical presentation of selected variables 
at the rotor inlet, rotor outlet, stator inlet and the stator 
outlet. Figures 7 through 20 provide examples of the plots 
avallable for display to the individual on request. The 
user 1S given the option of terminating the plotting sequence 
at any stage of the presentation through the use of inter- 
active prompts. 
The following is the subroutine's algorithm in out- 
line form: 
eelgorithm: 


Convert the appropriate variables to Real 4 for compati- 
brity with the library plotting package GRAFF. 


Determine the values of axial velocity, relative flow 
angles, total-to-total pressure ratio, and the adiabatic 
efficiency for the rotor inlet, dispiay if requested. 


Determine the values of axial velocity, and relative, 
absolute, and deviation angles for the rotor exit, display 
1=f requested. 


Determine the values of axial velocity, absolute flow 


amyres, and total-to-total pressure ratio for the stator 
inlet, display if requested. 


68 





Determine the values of axial velocity, and absolute and 
deviation angles for the stator exit, display if requested. 


19. Subroutine SHAPE 
Subroutine SHAPE calculates the eight nodal shape 
functions for a given point (&€,n). The equations for the 


nodal shape functions are: 


N(1) = (én + €7 + mn? + En + En*® = 10/4 
Meee =) 2° eon) 2 

N(3) = (-&n + a a5 n? + zen a ee - 1)/4 
N(4) = (1 - n* - &€ + En*)/2 

NOS) = (én + 6° + n° - E°n = En* = 1)/4 
cM Gen —-° + Ench/2 

NEG Joe (= oi + Ee + ne = aan _ ane - 1)/4 
N(@) = (t= n° + & = En = 1)/2 


The values of the shape functions are stored in the array 
Ppueeand are returned to the calling portion of the program. 
20. Subroutine JACOB 
Subroutine JACOB computes the partial derivatives 
of the shape functions with respect to § and n and computes 
the elements of the Jacobian matrix, [J], for a specific 
point (&€,n). The equations for the partial derivatives were 


obtained directly from the differentiation of the functions 


Shown in the description of subroutine SHAPE. The arrays D 
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and E store the values of the (dN/9&) and (3N/dn) respective- 


meee the Jacobian matrix is calculated by the following se- 


quence of equations: 


Joes, 


Jel 


Jul Ay 


wi; 


n 

1) = (92/26) at (ON. /98)z, 
el 

2) = (dr/96) =) (ON, /92)r, 
n 

= (ozo = } (aN. /an)z. 
n 

2) = (3r/3n) - (ON, /an)r, 


Arrays D and E and the Jacobian matrix are returned to the 


Saeting location in 
mi Subroutine 


Subroutine 


the program. 
ERR 


ERR1 is called by subroutine INPUT if the 


storage limitation for Real 8 variables has been exceeded. 


The subroutine displays the amount by which the limitation 


was exceeded and terminates the program's execution. The 


user response would 
Sible or reduce the 
Ze, Subroutbane 


Subroutine 


be to increase the value of LIMR if pos- 
Size of the mesh. 
ERR2 


ERR2 is called by subroutine INPUT if the 


storage limitation for Real 4 variables has been exceeded. 


The subroutine displays the amount by which the limitation 


was exceeded and terminates the program's execution. The 
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user response would be to increase the value of LIM4 if pos- 
Sible or reduce the size of the mesh. 
23. Subroutine ERR3 

Subroutine ERR3 is called by subroutine INPUT if the 
storage limitation for Integer 4 variables has been exceeded. 
The subroutine displays the amount by which the limitation 
was exceeded and terminates the program's execution. The 
user response would be to increase the value of LIMI if pos- 


Sible or reduce the size of the mesh. 
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Vico Ulkes eAND DiSeUSs TON 


A. PROGRAM VERIFICATION 

Four operating conditions of the NASA TASK-1 compressor 
were used to test the capabilities of the programs MESHGEN 
and TURBO. In all cases the 63 element, 222 node mesh with 
an under-relaxation factor of 0.24 as recommended by Hirsch 
and Warzee [Ref. 4] were used. Selected portions of the re- 
sults obtained are presented in Figs. 7 through 76. The 
points annotated as “observed values" were obtained from the 
material published in Refs. 23 and 24. The values attributed 
to Gavito were obtained from Ref. 8 and those attributed to 
Hirsch from Ref. 4. A discussion of the predictions for the 
various conditions are presented in the sections that follow. 

1. Test Case 1 

For test case 1 the operating point was defined as a 

rotor speed of 50% design speed and an inlet mass flow rate 
of 107.6 lbm/sec. The author was unable to locate this spe- 
cific operating point in Ref. 23 or 24 and must assume that 
in Ref. 4 the mass flow rate was modified to conform to the 
end-wall boundary layer scheme described in that reference. 
Therefore, it was necessary to use the observed values pub- 
lished by Hirsch and Warzee in Figs. 21 through 30. The 
relative differences found between the present predictions 


and the reported observations are given in Table l. 
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a. Comparison to the Work of Gavito 

Figures 21, 22, 23 and 24 compare the results 
obtained by Gavito with the predictions of the present pro- 
gram. A significant improvement has been obtained at all lo- 
cations. That this was possible was due in large measure to 
the solid foundation to the present work provided by Gavito's 
program and to the excellent documentation given in Ref. 8. 

b. Comparison to the Work of Hirsch and Warzee 

The predictions of the program TURBO compare 
quite faborably with those of Hirsch and Warzee. As shown 
mmeeigs. 25, 26, 29 and 30 the predictions of the two pro- 
grams have almost identical average relative errors for the 
velocity profiles. The predictions of Hirsch and Warzee 
tend to have better agreement in the rotor and stator tip 
regions, while the program TURBO has slightly better agree- 
ment near the hub. The program TURBO's predictions had a 
3.7% and 2.5% average error at the rotor inlet and exit re- 
spectively with a maximum error of 4,6% at the inlet and 7.0% 
at the outlet. The stator inlet velocity poredictions had an 
average relative error of 2.6% and maximum error of 5.8% and 
the outlet predictions had a 2.0% average error with a maxi- 
mum error of 6.0%. The prediction of both programs for the 
velocity profiles show excellent agreement with the observed 
values. 

Hirsch and Warzee's program consistently produced 


flow angle predictions with closer agreement to the observed 
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values for the published rotor outlet angles. Though the two 
programs had the same average errors of 3.3° for the relative 
flow angles and 4.8° for the absolute flow angles, Hirsch and 
Warzee's program provided better qualitative distributions. 
The difference is clearly shown in Figs. 27 and 28. 

Hirsch and Warzee did not publish predictions of 
total pressure ratios or adiabatic efficiencies, so the pre- 
dictions made by TURBO for these parameters were not pre- 
sented. It is assumed that no significant differences could 
have occurred because of the similarity of the results for 
the velocity profile and flow angles discussed earlier. 

2. Test Case 2 
Because of the apparent modification in the mass flow 
rate which was assumed in test case l and the lack of compara- 
tive data for all guantities predicted by the program TURBO, 
another operating point at 50% design speed was compared. 
The operating condition for case 2 was defined as a mass flow 
rate of 114.7 lbm/sec at a speed egual to 50% of design, 
woeeh corresponded to reading 38 of Ref. 24. The results for 
case 2 are presented in Figs. 31 through 44, and a summary 
of the relative differences between the predictions and ob- 
servations is contained in Table 2. 
Bee lest Case 3 
Test case 3 corresponded to reading 45 of Ref. 24, 
which was defined as a speed equal to 70% design and a mass 


Flow rate of 151.55 lbm/sec. When reviewing the results 
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presented for this case, one should note the decline in the 
agreement between the program's predictions and the observed 
values. It is the author's opinion that the degradation is 
primarily caused by two factors. The first is the formula- 
tion used to compute the meridional velocity change across 
the rotor and the other is the application of a single block- 
age factor to all nodes to account for the end-wall boundary 
layers. Both factors are much more significant at 70% design 
speed than they were at 50%. At 70% design speed the rotor 
tip relative Mach number is about 0.94. This would require 
the program to account for transonic effects at the tip. 
Second, the mass flow and absolute Mach number of the flow in 
all regions iS Significantly higher at 70% design speed. 
Therefore, it is unlikely that a single blockage factor will 
work satisfactorily in all regions of the machine. It is 
hoped that both areas will be addressed in any future work 

on the program. 

The results for test case 3 are presented in Figs. 45 
through 58 with the corresponding differences between the 
predictions and observations summarized in Table 3. 

4. Test Case 4 

The operating condition of test case 4 was at 80% de- 
Sign speed point with a mass flow rate of 174.54 lbm/sec, 
SemLesponding to reading 50 of Ref. 24. Figures 59 through 
72 and Table 4 present the results of the program's predic- 


tions and comparisons to the observed values. As expected, 
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and for similar reasons to those cited in test case 3, the 
agreement between the program's predictions and observed 
values is significantly poorer than any of the three previous 


Cases. 


Beef OLNTS OF INTEREST 

The results produced by the program are highly dependent 
on the value of the general blockage factor used to account 
for the end-wall boundary layers. This factor influences 
both the quality, in terms of agreement with observations, 
and the stability of the solution. Figures 73 through 76 
show the axial velocity distributions for the rotor and sta- 
tor for test case 4 with a blockage factor of 9% instead of 
the 6% factor used to obtain the results shown in Figs. 59 
memeeugh (2. A comparison of corresponding velocity profiles 
clearly demonstrates the factor's pronounced influence on 
the program's solution. The general blockage factor also has 
a strong influence on the program's convergence rate. In 
some cases the factor can cause the program to become oscil- 
latory or even divergent. 

For the low speed cases of 50% and 70% design, additional 
blockage factors had to be applied to the duct element be- 
tween the rotor and stator tips to obtain accurate results. 
As shown in the program listing for subroutine INPUT of pro- 
gram TURBO, the factors used for 50% design speed were much 


higher than the factors used for 70% design speed and that 
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no additional factors were used for 80% design speed. A 
meebal method of calculating the end-wall blockage factor 
must be incorporated if the program 1s to become independent 
Seeinputs other than physical constants. 

The deterioration of the program's predictions with in- 
creasing Mach number and its failure to run for test cases 
with strong supersonic relative velocities at the tip demon- 
Strate the need to provide the program with a method of 
handling supersonic relative velocities. Hirsch and Warzee 
[Ref. 15] showed a method of extending the radial equilibrium 
formulation used in the present code to supersonic flow. 
They presented comparisons of predictions obtained by this 
method to observations of the NASA TASK-1 transonic compres- 
sor at 100% design speed. The results were impressive and 
clearly showed that the method is valid for relative velo- 
cities in excess of Mach 1.4. It is strongly recommended 
that the first effort at improving the code be an effort to 


modify TURBO to include the method shown in Ref. 15. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A computer program based on the finite element technique 
has been developed and has been verified satisfactorily for 
computing flows through subsonic axial flow compressor 
Stages. Minor modifications have been suggested to allow 
transonic stages to be calculated. 

The code was written in such a way that it could be 
readily adapted to compute either turbines or compressors 
with multiple stages. Before such extensions are attempted 
however, the following specific recommendations are made to 


improve the present compressor code. 


A. PROGRAM MESHGEN 


., Incorporate some of the two-dimensional techniques of 
Adamek [Ref. 25] to improve the efficiency of the 
code. 

7 Review the code to find improvements in storage allo- 


cations and computational efficiencies. 

oe Modify the program to track the first and last nodes 
of the rotor and stator as subscripted variables so 
that the program can be used to generate the appro- 
priate mesh parameters for a multi-stage machine. 


4. Convert subroutine MPLOT to the DISPLA system. 


B. PROGRAM TURBO 


is Incorporate a method for the global calculation of the 
blockage and losses created by the end-wall boundary 
layers. 
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his 


Convert the storage of [K] and the solution technique 
for the program to at least a symmetric banded scheme 
or if possible to a skyline equivalent scheme. 


Test the program on a variety of machines and operating 
conditions. 


Obtain expressions that approximate the NASA correla- 
tion curves for 65-series blading. 


Incorporate methods for the prediction of stall/surge. 


Take advantage of the modular form of the program and 
include a variety of correlation techniques as a user 
selected option. 


Review the program for improved storage and computa- 
tional techniques. Specifically, determine ways to 
take fuller advantage of the dynamic dimensioning 
scheme [Ref. 16] used by the program. 


Convert the rotor inlet calculations to allow the 
value of 81; to have nonzero values so that the program 
can be extended to multi-stage analysis. 


Modify the use of the values of the beginning nodes of 
the rotor and stator to subscripted variables so that 
the analysis can be extended to multi-stage machines. 


Convert subroutine MPLOT to the DISPLA system. 
Develop iterative schemes for calculating the flow 
angle, the velocity distribution changes and the 


thermodynamic property changes across a turbine rotor 
and stator for inclusion in subroutines ROTO and STAT. 
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Soneect ly tly Relationships for Figure Three 


Element Local Node Global 
Number Number Number 
if 1 iP? 
1 2 8 
1 3 i 
1 4 2 
I S 5 
il 6 9 
i 7 14 
1 8 ils! 
2 i 14 
2 2 9 
2 Ss} 2 
2 + 4 
2 2) 5 
2 6 10 
2 7 6 
2 8 ie) 
3 i 16 
3 2 ne) 
5 3 5 
3 4 6 
3 5 7 
3 6 sles 
5 i is 
3 8 17 
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TABLE 2 


Someasi son Of Program Predictions with NASA Task-1 
Compressor Measurements at 30% Design Speed 


Average Maximum 
Difference Difference 

Recor Inlet 

Axial Velocity 4.6% Laas 

Relative Angles Po? Brae. 

Total Pressure Ratio Oras 1.4% 

Efficiencies 4.3% 11.4% 
Rotor Outlet 

Axial Velocity Se % 6.73% 

Relative Angles eG * A 

Absolute Angles oa 2s 

Deviation Angles A as 30 
Stator Inlet 

Axial Velocity 3.4% 6.8% 

Absolute Angles 25° a ad 

Total Pressure Ratio 0.3% Le 2s 
Beacor Outlet 

Axial Velocity eas Dac 

Absolute Angles G6” ioe 

Deviation Angles 1.4° ZnO 
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TABLE 3 


Compamison of Program Predictions with NASA Task-l 
SOnpoecssor Sicasuremenes at 70% Deston Speed 


Average Maximum 
Difference Difference 

Rotor Inlet 

Axial Velocity 4.8% 7.4% 

Relative Angles i se 

Total Pressure Ratio Leas 2202 

Efficiencies Be Cae 4.83 
Rotor Outlet 

Axial Velocity 4.1% 7.4% 

Relative Angles 2a Siar 

Absolute Angles So oy ok 

Deviation Angles EL Bley e oT ae 
Seator Inlet 

Axial Velocity 4.0% 8.8% 

Absolute Angles oe is 

Total Pressure Ratio 0.5% Zaaee 
Stator Outlet 

Axial Velocity 3.6% 7.4% 

Absolute Angles 0.6° oN 

Deviation Angles eo 234% 
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TABLE 4 


Comparison of Program Predictions with NASA Task-l 


Compressor Measurements at 80% Design Speed 


Average Maximum 
Difference Difference 

Reeor Inlet 

Axial Velocity aos 11.4% 

Relative Angles ee Se 

Total Pressure Ratio 1.5% r 

Efficiencies 3.1% % 
Rotor Outlet 

Axial Velocity GOs 8.9% 

Relative Angles 20% 5 ee a 

Absolute Angles Gio Sie) 

Deviation Angles Zea? Aa? 
Stator Inlet 

Axial Velocity 629% | Bs eh 

Absolute Ancles Deo Ti Se 

Total Pressure Ratio 0.9% 3.1% 
Stacor Outlet 

Axial Velocity 8.0% 2020 

Absolute Angles Oars = 1 ee 

Deviation Angles eS * Zu 


86 





Pigure 1, 





Sl and S2 Stream Surfaces 


oT 





juoWeT| epON-3ybt_ ue FO sanzeTOUSUION “C S105 tat 


L “ © 
\ 7 \ 
\ | \ 


| 


c UOI}EYS 


88 





YSOW JusUeTY soIyL_ Pe Jo otTduexg 


eo anbT 4 





62 





aueTd a’z ayq oF U’3Q ZU WoT SAtyusuStae lea hutddep 


Z 


JUuawa[ta Way, TFAAND 


| L 


quausta yuasied 





7 \ buyddea 


"p oanbtd 


90 








on 


Mesh Generation 


PAguve 32 





MOTA-UuhNoOAUL TRUOTPTASW SUF JO SuOoTReTNOTeD sy 
UT PeSQ USOEW SPON ZZZ 9A FO 2OTd paqVerSusy NAOHSAW *9 oinbt gy 


(SoyoUuT) BsoUe_ISTG TeTXY 


OU Ol obeepeC rect 0Ce Uc oc ve ce 0c Bt ol rl cl Ooo 6 9 Gat ca 0 
—t71-+-4+—1 0 





>—-9—--—_0--—_—-9--- 0 —_9 —0 
| 





Be 
fy 
0f OQ 
}- 
#21 © 
a 
vl U 
oh 
et py 
S| 
0¢ 0 
0) 
ce 
HT 
pe 3 
92 4 
8c Ww 
0€ 
ce 
be 
9€ 
Be 


© 
HSI INIA IIULIvUVdOST OOM 11973—— Op 


eZ 





Coleco coc Vet acu 82 
Sn oa Sea ae ae ee 


os 
7 ee 
pee 


: 
* 
Docs 





° 
t 
‘ ® 
e ® 
‘ 
e ° 
e . 
. ® 
° 
s e 
. » 
reg *& of 
. * 
’ 
° . 
e . 
R . 
e . 
* . 
‘ s 





ad 


RA OI a ela balan Stata eee cecee cece sera cecececerecses sores c-cccccsccs-cecseenees corey © 
° 


° 
; 
° 
i 
° 


* 
e 


oe Soe emt age tase cesses see oa cowatcerone-eobe OR OCCT ESS 0e CoHe C8eS ee et eeeee He 
ry , . 


ry 
‘ 
cy 
e 
. 
* 
e 
e 
e 
® 
® 
° 
e 
* 
° 
s 
e 
*> 
* 
e 
e 
e 
e 
e 
» 
e 
° 
* 
° 
* 
. 
° 


© ce- emecee: cecGeececcccsccoen chon sceecccceteesd 


20 ome $0 Same O07 oO SSE EET CUT WEES gE S SEE CEES ee cece 


eeeeeenes 294° 


; 
i 


poods ubtsaq 


(990S/ 


iD Q4 
4.) 






eee eee Sree 


. 
° 
+ 
° 
’ 
° 
. 
e 
Tree rece cece Moca en cecccncs 
iJ 
° 
. 
e 
e 
. 
° 
e 
° 
° 
oo > 
. 
° 
° 
ro 
° 
. 
: 


eeceet ares 
= Se CeO ee eH e-Peosseseeseeses BOs cese.2 
ry 


’ 
e@e +8008 ©: 90 -@: e+: 8 


e 
. 
. 
s 
° 
° 
° 
. 
° 
oan ee toc Seas oP oosesevece 
. 
. 
e 
e 
® 
ry 
. 


Phe merece eee © ore © Ser ce cee re cece gaeceesses os 
‘ * 


’ 
«Ce wecc cet ecen-coceeens- caress ce ces os 
: 
‘ 
oe ee Sere ee eee Serer erie er rey 


eceecece 
ee ee ee ee re ere Teter ee 
i) . 


svt cemmeccvsterbooscummewcsce “wt ee- eee as 


’ 
’ . 
Pee eee eee ee ee 

. . . 
. 
e 
e 
e 

OO + 08 +s C8 OF Oe He FSS S: OFS wee 

e 
* 

. ° 

=Peee cogces stors+ -Gesessese®-eseeee HF eeven 

. 
° 
ty 
° 
. 

Pere COMO TOSSS SOC H ASSESS ESe © 

. 
° 
J 
e. 
e 


ee ce P- ce. eee oproee- 


+ OOO 0e- | eee Seacces- ©: ww eecere-res 
. . 


eoecege 
eT eee Tr a es ee Pe ee eee 
. . ’ 
. . 


eee see of wee 


+P oreene 
2 eer accoe 
va eeeeseserooe 


wh OB FOO So SEO OF OCF ORSTES OFS COFFS MH SON TET SS SS SBSOREEBEt SC. 288. 2-6 -: 
ry 
e- 


sbae-eeet eer e G00 Peasseces 
. 
. 


* 
—o- 
. 
° 
. 
° 
’ 
° 
° 
° 
. 


ee. Coe 
sererces 





: 
,eeee 880004 le eee secee 
aRgiiceen 
i 


t: 
+ 
+ 
‘t 


’ 
* 
° 
r 
{ 


‘peeseescsce ames 


. 
+ OP cOCcemmee ee eer ose 


’ 
° 
Sco eeeewecccecategen ses 
. 
‘ 
’ 


conse SOO © 'eeere 
t 


’ 
ep OOF SCOOPS CS SOOT SHRSE HS PP co OmPeses SHOFSES - FeEnGES 
° . . 


. 
° 





° 
° 
° 
° 
. 
° 
. 
® 
ee Pr 
. 
‘ 
. 
° 


t 


po--erecemmene- ope er Onespescerbocanccscccoree-b - 


0 
° 
id 
ry 
ot eeee coewes- f 
. 
° 


§ cc eens coc emamecfe os: wae - wens copeters 


. bac coene oe om 


‘ 5 
oe cee Cn ee ot ag 80 cc cme webs Ge Ogee CONS HEED OT HOSES HOES EE TORO EEe BS: 


woweames = cos etebe-coeeseese o* beccecewne cee ss cheo--cce- 
cow cle cewce se ww ecegacescec ere: Ba cetowes- oer 


cotwecnecenne cpesenmeccereesne pare 
eee eee ee rT 


So cam +m te cece eg cesses 


come eneogor™ seehoevece eo Come sec apesure ceive eece t eacecr-seteacr 


OSS 04 2 8OO* OS OO eh 8 OES OF OFFS WHOS CFB OO 


Gann ceaupaneancbeccmSsvce-0ac0®-<o=teeene ee 


Pores mamas cee eh- an -ccncccon. > 
: 


die lak PigcSnOAD SOS Ocoee 


qe AQTOOTAA TeRTXW OTUI 
1OVOY SYR JO AOTd STO XTUOFYSL poyesrsusy OGUNL WV 


TOOTSA TeX 


en eB oe 


recOrKy e Olesen @ -¢ 
ery ¢ 8-8 06 See ie 2 
2 Geeta 
ee. ot fe 
» . e e * 
s 8 e ° ‘ 
bey at pee te nance noel ae 
6 en Ge-ere free ot ce Ce ae 2 
: : : ° > : 
es gee te! Se aha 
ites te 
° ® ® . e ‘ 
se mateect sae 2 
yet es 
ti oe a oe oa 
eee 
° : ‘ i : 
sive guaisl cone 


ove cahcovecescocecce pee: eeereece Snes Pa 


ote a? a nen oo 
: 
; : 
: H 
enone eete amo e Goce 
i t 

: 






eotsece Oe Ie 
‘ 
‘ * 

Per ee He mew of e 90 os case ec cee -Po suse ecerssace- fo > 





s s 
Go one og- © Coches ce-qre-ce 


sal 
| 


2 e+eeen eres 
eee cet etco= 
99 © Oh oweee Hee 





ee Te 


. 
° 
tereecenher * 


) 
° 
’ 
. 
. 
. 
e 


‘ : . . 
FO CeO OO SOB Od 6OF SF Fe CORRE) OFs EMSS EOFS BOs See easeae: 2: 


. : * 
Se ee ee a ee ee 
. 


= POCO: OF =o O08 5998 0SHS OEE apbeeseoemeseeccem:: ++: 


. ‘ . 
Ce ee ee ee ee a ee ee eee 
. . . 
. 
° 
’ 


° 
i : 
ceece ce eee ah -o 8 eed eens on ad cumece ce sec cecrbes 
‘ 


mde ce cc cmc cw es ony oc mac scee ess compose: -o8se-00 


. 
’ . 4 
POEs BLO wh & OF CERES OF OF 6 © OOP SOOM S FSD Ss OOS GED OS © OS HSS SOE OF EB G0 OS FO SD 6 FOSS 6A SS OSS GOSH 2 G0HS BH Os Seoe-= 


. 
© cee mweee gs ah eeemmeecce® cheseen® 





BCE OE Ores 
° . . . 
ee cwecsecsestes COCO CF COFCO CEES + Bs SESS 1 OFS OOF = COBDS OF SEHOCOSMMOSe SSOP ETROL ES Masses ees sersete mbes -s- ++ 


; F 
haghh ooo ete 0 00 oe 00 0 a oes ene = BB Oe COTO ESOS OOOOH OT ES CHES SOS ER OT OES: TOM OSes COCO DOSES E> OBS es Ore HOS OEE SESS ODOT ese eeemme « 


$ 


ey See 
— ; 
o- Peete re ve 


° 


oo 


", vanbty 


Rerves 
> 
~~ 
—t 
med 


™s 
Tpea 


Ta Te 


aoue ys 


(seyout) 






ie eet 
Sonar ca 
oe 
i-4 4 461 
ee 
ie 
rot GPCI 
oe ye 
i> he sane 
eae 
ini... £8 I 
| eee eee 
j i 6 
ne 
wid Ley 


o3 





69 89 


SOTHuUY MOTY SATRIeETSY JYSeTUI AORFOY 3y 


JO 490Td 8T9 XTUOFASL PpoxFetouSH OGUNL V 


"g oanbty 


(Ssortbaep) s[Tbuy MOTA SATRIETOY FOTUT 103904 


£9 69S 





Boye ude Co (8a eo oS. BS SG 8G “SC 7c Fe 87S) 1S. CC ar 











: eeete ; Sane Or ee + ence aes [ . 4 . ee, Ce ee fo d | 
“i onj--- ere rower. , t 7 q a0- t ? if i 
t : 
: 
cas t- : = t i ’ BD ccceere: t Secmnesice Pooseccee. > C | 
. ; ; 
cae ae 
= “e- —eoe a *e t ] . he es oo et: o-4 o 7 smetneces € if 
: as oe en $-- oh e0g-e- co oopee oe See eaace Moor Se et 05 Gee -odadsode JecCuOCGOOr SaconoFod ) mule O) GSU ae lnnder? Se Soenasoerl an pesmoce ya eoeeG. O aeseesaee pl 
: : : ; ; : : 
: : : : ; : 
ee ak ea 7 
: : : t : i : 
bd ommesoes Qaveces can ae @enaces ~- eoee atee Be + ceeee 2 eea-eetes hro-ee . : - 
i : : i St 
: : i 
: : : : : : : : : 
HY e ° ° e . 5 5 
: : i : : : : : : 
es Pe einm a psecsasess Bietrnotens $ econ ey eerate @ecceat oars . , . > . = : M : : : 
: 7 3 : 4 : * : eccene t eoeaet-ae x eeccctoce ge alae Lee edie ae MLN ok! ba oo * reer 4 - cud ee oom ° = . of ececce: 9 T 
a See Se 2 ee po Ul a : UE 
‘Te ek ee ee ee ie oe ee 
i : ; ; : : : : : : : : : : : 3 
ee eae eee ee. ee: oa 
- OO PEO lOc oP e- Cees es med cnet es Po Be ere ne te teccece trees aoa wetech eves s cong saeseeesahs saeseor ss eceneees freceemensepteasenereapersmnscee gee ances fee see dee be! wargeee-cee see] . l af 
H C s ' ‘ e ° e = ° : : : : 5 ° O ° 
; ; : P ; ! F : : : : : : i 
‘ : I ; : i : : : : ‘ : i ! 
° ‘ e é HM 5 5 : : : ' 
sats SE ee en Oe ee ee |e a. or a a: 
Ee en Caos I ee eee A ee ee ee ides eae | 
: : Ween terme ore eS ha ooo} © sweeeaes Fonnerem- fortes erceecneemeneperereccee fs mrmeacenfseceeseeepemtesanee fees sear eGace casas hea tee fr ocwcese a eeeae g | 
. : i : j : H ' ; ; i t : : : : t i 
4 : : : i H ; : : ; t ; : ; : : H [4 
: : : i i : } i i : : 3 t ; : ' ; 
: : : : : | : t : : : i : : ; é sitios : 
s i : : H . ° ‘ © * ° e . ° e 
2 aoabeut ereces ee etaunce OR OICOCHOOC SCRE ree ee ee ee ee ee on de WE Qo10y Jy $2 PNY ROIS BPI VIFAA—_ 6 T 





soueqstq [Tetpey 


(soyoutT) 


94 





SREDEYM PeTuy SA«t07 Oy “oUt 7O-OL4ey -oansseig 
TeIOL S42 JO 3OTd BIO XTUORYSeL, pexyerzsusy oGuNnL Y °6 esANHTY 


OT}ZeY SANSSSeAg TeROT, 10704 
Gola lL ewecerme GG. “Sel. 2cl : Oren aac oe BOTT por: 






Nan ses qratce- Greece eee onegres ee eee ere Gute F teee Preseene Ecotec: Oo ee see oO seer Betctors o* fe ee a e ot coce * . ¢ Py aT te { ’ 0 | 
Rice cn0G. rncne Fann ncheeee chactaachinsc-espboconag soe ce Peet oe 2G ome Poco teen e eaten ona Sores itee ee Ses eo Ce Coc ceo cers eens t to e Soe eee E T 1 
ee a) ee Od ieee Si ae : a — ae 
te eee) ie cee ener Oe COD ec Se tte ae Oe re ee ee ee ; 
ee he Pett ey a Te ie Soc Ne ae > ey. 
~ see§eaccs ‘ avacaguecessy ee >: ° pavers acne ie ee t- oe cote ate spaces men setemmecneeg sateeg cones fee ‘ soot t- * t “4: < seo @ 4 | ’ 4 é | 
yur, Ste? <e re ois Boe Sa Se ee ete : 
: : x ° ° A ° : : . : : : ° : : : : : 
So pe ee So me os che OF ; ey ee ety : 
ee Oo eee eGo) eee es) an ec Ne ie ge ac ee Sc a a SE S&S 
we Been ter gsaeuanseeeser| weep at car teenage er aih esene de comeege moot soneefeosoeeh P rere thectes Fema eens fee : meee Groce: @ t aa Ponce 5 aoe et: ' eo. ¢ t ' e T 
: ° : ; ; : : : : : : > : : : : : : : . : ‘ : : : : : : 
se Gs) cepa g Bete Begs Sta et rit tee i oe: ee a i 
St ceca a ait foe re Ont see semen sh Sy simmer eae Net et a ee SY ae 
® : : ° 2 : H c ° ‘ . : ° : ° 5 . 2 : : 5 : H : : : : 
Ss Satan ares et Wis Gt | ict ic” sae et ea ae et ee ee 
bs Prony ie ets eon f alld oad ' Po-sfeee. ee poe Wgahr ores euarecrcta we seheese ‘ ee oer igh - bla Jere ie ere oot ' gan *4 Shae . *t . t t t vi 
a es ee: ee ese, igo? at Me ee ce, cee ee tee ee OM et Re Or ee 
EE SOOM Sae Cas eee ee oie No: ik MMM? VUES Wiha VN tie SE NG aht= oe) ae 
ee a ee ng ee ae ee a 
ON Cine ita ste re edt: pO arn aed thy eae OME Rene! Ct MSU A ciate Rl Sete eee. Cas cee ee 
Uae abate Goals co eet came. Pec cobe wt Pownce Powe ees omnes Fo comme ho cease Porves tence oceae & patted w eterg © ae | retseprotes r ccs Boats Si ae oh Se fees ati : - t on | 
ee ee ei os a ee mee ee ee es a eee ee ee 
Sue ee rte te ee Ot eos eee CO Eur st ety ON So ep gee cea oe 
Re en NR es “ai 2. SG PE oe Rs cate St Gg ot este el see) oye 
: : : : H ‘ ° : ° H t ° : : t : ‘ : : ° : H ° ° H ° ° : t : : 
5 e ° ° a ry a e ° 6 rf’ e . e e ’ e ® ° ° . ® bd . : . ry r e 
° ° : H ° : H ° H : ° ° . H $ : H 3 t : j : ; : : : : * : : 
s 4 ¢ ry : : ° rf ° rt e ° . * e e ’ ® ® . : 3 H : ' ® K H 
SOM Re Ne eee hee a Codes Cte er) okies camp epee la heer ct CaN hokad soe Pa cy et ee 
ale aie ees WA Na a ada rca i aa, doe ae gr aad A Pl aaa | swencpoe-tecfoce eee Serres Ct. a ad t woe ao te -9 1 
: A : : ° ° ° : : : H : : : ; : : ; ° : ° ° ; . ° ‘ a : ® ° 
SA a Ie ke a a a de ee ee ecw 
Sau are eect We ee pe te ee ey ee a eee oe ee a | 
e e ° ' s . t ° s e . . ° . ° e ° e e e . * r rT . . R . iF 
i ° ‘ : t : & : 4 : é H : : 4 ‘ H : : ° 4 : : : : : : : : : 
: 2 4 : $ : : : H : : i : : : : : ; H : ‘ : : : a i : ; ; 
Re yet oa See os Se ete ee ae a ae. ee Oe ee On eee pea. purists. bet) 5-5 eee 
Shs seer eect eat gage test eeoere a each ain aaa deel adcrach aio Beat aw eee paid ie ces ase : ° se i - l T 
ee et ones rs Se ek O OOe ee eee eee Py Ce aie” Fee ae: 
ee eee Por eS me ome i UME lS Be) gn: ue 
° ® e ® e | e ° e : t e ® ° , ° 3 e . . * . : H ° ] H rf : 
ee ee te i es a a Oe iene ee ; ot > 
Si ee aes es eel. fe, ee LE or ere ee ee | 
sa ag ee ee a tt cate aes el es 
ies ehh ah rao Sas ay Si Sin lc Gn hens ics nats ey Sangetnh arc gees ee azole ge SE cis Gsm § aaa bigs corsa tae Q 1 
fe re ee ee ee eee ee ee 
3 3 5 ? 5 : J 4 ? t ' * ry ® t : e ; ry e . ‘ : 8 ‘ $ : 
: ae eee eee ae ee ee Se ee eee ee |e ee eee ele we 
Se ee ee ee ea eae 
H 4 H 3 : 4 : K . s 4 p r : e t : ‘ : e ‘ ° Pe i i i : : ‘ 
-_ loot tot oto beet db dt gd tt OGM. (DPINN aN NSSsudnmei gia 
Se mteeeFemena+ mecie sieSerlece-#-nsan ein nlcin dM oveuseh sense Su lemrs@osensabertcce-seereschsecctatces cocmre cecum cere mdse cee ean bene ene eee n eens oes See ee ee ee ee ee eee os es a O1 Ore 


s0ueqIsTqd Tetpey 
5 


(seyout) 





SNTpeY FOTUL SA AORZOY 9yA Jo AduSto ys, 
OT}FeGeTPY SUF FO 2OTd BIO XTUORYSL paezyersusy oMING VY “OT eanhty 


Aouototssqa OTFeqQeTpy 10704 
(6 96 (50 SU ome ccs ye oe Sp CCte Bee Ge Sere RRs Sr ce | 








soueysTtqd Tetpey 






eee ea a Re yy eee ee all 
\ —_— 
me 
tae re nae ae ROPE eA eee : SI 
i 2 
Q 
Ss 
i : i 4) 
. ae @ wee me ob ee eee | Toad ox —-o= "* Q] M 
: : ° : — 
; : : ; : 
3 : ; ; 
: : : i : ; ; 
em ee | ad ce San per a 1 
! : : 3 . , : H 3 
: 
. : 5 : 3 e ; Ld 3 
: 
= “a ase? nm connetnennssenesem§-renasenes can Ge -ance—a aA oe nema Pree r ewer enen ee hs ae aaa a serene cee age mee 8 | 
: i : : : i ; ; : : 
i i ; : | i : 
? 3 ; i : 
| 


asians Fe ee ee AW WOLOE Vv AW3T9T 993 JLive vi Cy a — 


me ee +e tame Bene ceee cme et ate sete ant em ommners eo: = 


96 





ATOOTSA TetTXY 3eTINO AORVOY 3yA 
JO 30Td BQ XTUORYeL pazerouay oMUNL WY “TT sanbTtg 


(9o9s/ajy) AQTOOTAA Tetxyv 3FaT3AnoO 10304 
Bvt etUS soe) Gur OVE Cll  seUe che CEE SOC ErSeh OCG. oce 2Cl Oce Sle S-6 rifle 
' CT 


4 
















































+ ? + co 1 s on ' ] (i 
: : : : 
—- 
ae oe 
: : Sanne Pete Ere pint etee wrt ; : : ae on) a 
bah eras sh glerigis oh gas t oe isa ia 2g —@¢ a= lak hares ces ica ‘ ‘ ¢ 1 oe + . eef--- eee t t oe rang ri T 
: : Pan: : : a oo ria eae as aes Goyer: : ue oe Bees 
— ae ee te Ste eee ee ae eee a: 
: : Ee te Otc tye Ba Gh Cee oma Se Guwel stk ae derey ks rg ade Phe : 
es Sears ts. de yt Des ea eS SOR del et heme ety Rita k Ged Ser: eas Ne 
: - 3 eek ee ae Seite Giect «colt Giri mmun es et s dueien ts TS oi : 
tie a eo ne ae eee ee ae St er 
ae —fe Ma ee f e- t ee rod => of * aah rics ipa ane A rt le ah ae aL ee pecs neg a os dh ? wee ! ee t rs i sane t 1 
es ee ee a ee een Ss Mb age ee Nee ee oR Ce a ts ee ee oe 
e » ° ‘ r) 4 e ° * ry ‘ ry t e : e ‘ ° ¢ ry ® ‘ e r . = H ; 
Sa ie ane Se es aioe te teeenti suet mee nas sides aioe a: oie Ler ee ae 
tr ee ae ae Are isehuionit ed Cae pa er Ss ee See a: er i ee eee et oe 
Be ee ma eon oe toe ie ete, te Ok ae ee ee ee ae ee ee ee : ae 
i Se Re cae: ae as eae ae ae ene Sy CN oe Ca nen come ee 
paste eits es > lion. 4S: ; 3 ae ec 0 ee ee ee ee eh ae 
de ERAS patent SO: See eee oe ae ns eee ee ee ee ee ite ae oat tn een pee 
ae gee ase cee SE re ee ee bai aad -} mg tiate  Beasabtin dd wm eat moo Go- oh 8 obs aw og: ¢ nae Hoek fe mcePoce t=} 4 oe oe es Sa 2 vl 

ae ie ae ae ee ee ae oe ee RS ak wes 
ee ee er a ee ee : : a eres ye mre aro copes awe 
oa eH Ate cee tn eee Same. ae eee ee er ee a bof : 

e . e ‘ = a 5 2 ; $ : = . ¢ . e ° ° s ° . e . . s ¢ . 
- 3 ay ae ee ee 3 aes : ee ee eo ee ee eer ee peer ie 
ae aie ae amet tg fe eee ee ee ee : beri ois ws om 
ea : ee oe caer eae nt : ames eee : A ee eee eee eee. Ol 
hea ane MM arte cht ace ras) > : ne eeaerat cht Guth os bot ct sage ee tac ar@he we teean ae, 
Spas ae ran ed ee en ee er oe ee has Wee eA errant PGs ecegniG Wei iia omatee gc. At: Sten Piet 

ae eNrers ec0 aber eb anas eto ST , mies obese ae batemes + Qeewep agce Poece ccc. be meng. 2 Geo 4-- rere ge p- | es OR ao a t c C- Gene § -G { 
3 A : : . : $ : s : 1 : : . : . ’ H : H 3 : H : H ° : ° : : ; : : : G 
ip res nr ne egies © Pete (ete wth STN Lore, Oa eee ie ALE NPee ig eO) Pie oe ghee Sah eek «Us se 
OS UR oye = gs 0 pike Wir aoa eet ecto a an ee eee = SCY Ok es eee Ser eg ea ee Li toast 
° ry ° rf 3 ° : : 5 ‘ . & | S be : « e ° rf a ry { e . ® Py r ‘ ry a Py ry t ® 
' : 4 Py ® a 5 : t : 7 4 A $ $ C 2 * 8 ‘ ¢ e e ° . 3 8 « » * « . a e ¢ 
oe ees Reet rt Weer eer em Ore eigen gee eh len et aes a ee Re oe NAA es) oe 
a e Py e e e . 5 7 . . Co) 4 p g . a ¢ ‘ $ . ° 4 « * a PY 
a aes Ge 1S Ls Va ters. ae Fe H ge as ee ep ie ee, Mee eet Se ee ge ce pes Ra BE a ps 
Se gues Merry reser "ice 84g. oe topes. fermi tirrceis 0 eee AR Gee Mee eee Gey “ig 0 Ns GEN Beata ye en ee oes 2 tee oa 
a ee a REE ns ae Dates Cn pe ews lar me cee ay GM eee a ee oe a Cae Ne ee ee ee ee 
: } Day : : : t ? t ee : ws eae ams Ga aaa es eae ae ate Paes eae a ge i ales eae “fap abel ecole ape ala eam -ne ea 9 | 
ee Pais ie ieee tn Rae: eet timer Mime ccest ea er tk ples Peo RS ae SR oe 2k ee 
1 4 : : B : : H ® : : : : : i : : ; 4 } : i : : : : : ° : : : : : 
1 : ‘ ® . ® H ; $ : q t j 3 . . } : ‘ ® ‘ H 1 2 * ‘ t a ‘ 8 ° ‘ Py ® ® a ® 
Ae ait fe eet ne ts Ot ete fs ae mee ee So ee ee eer aes epee et ee eee en 
He og ie i ee be, on a ected Ad erie dase et Joe tig eases OO Ee eS Sa 
Coos eee 0 OM i Pt ie em kar MEN MGs ir MSU Oh Se tae oe 2: ts nee ae ae me 
ee ae ee eee er eS et | ea ae ae eg aes Pe one ee Ly eee eee oe ee ee ee 
° a e * ry . ° Py ¢ : i . : § . j G is $ : t a a ‘ ‘ ® ° ® 4 r) t ’ + 
Sa ee ee a me ee ee ey ee ie ee ee ee ae 
thane = a Scechccsc@ sas Oine8 s ~ S ae Op ae e ieee ae ry ° a ry a ° ® : : : ; Hy : ry : e : 8 ry 3 
: oe LG 3 4 : r oo $ 4 1 By t 5 mnnePian Soines 2 FF : a-ak ma ars ae bee OS t= a ae oot - > ea ct aed ss aa 
ee ee ee re Re ee a mee a me ie, ec 
$ ® 8 ° ¢ ° : . C ® I } « : r) ° : ° : t ° : ° ° : : ; ° ° . : : : ‘ ? ° 
4 a ‘ ° a ° ° t H ‘ i ° a ° e ¢ ° a 1 e Py 4 } ° ° ‘ r ° e ° e 
ae ee ee ee lo one a ere ee rs ne a a er or on i aie es rast 
ee eee ee rte eee ee i a ae ee a haere eee ee See eT ae oF ee 
a $ a e 4 ’ H : t P ° ? ® ® & a . r a a a : 8 : 4 rc ‘ ‘ ° H H i : ‘ H 
: : : : ‘ : : : ' : ; : : : : : - es ; H : : H : : ' : 7 
i ee eet re a re ee ee eS eee ee ee Ct ce) ae me": 
He ety 5 eee : aa ne ae ee I agate fr 4 ; Wixve: 

Se ee Ae Sh pi ede Rl ee ee AM UROL et ee io bf Eatin dorde a dso q ALTIONIN Wixv = RI 


SoueASTq Tetpey 
97 


(sayouT) 








SoTbuW MOTA SATIETOY 3FETINO AOROY Sy 


JO 20Td BI9 XTUORYSL pexerzouey OGUuUNL VY “cI eanbtTy 


(soartbap) aTbuy MOT A SATRIeETOY FOTINO A004 


|} $$$ $$ $ $$$ 





Pes commeac® 


se 
ee an 
a 
= ry ? - + ae S48 6 «ee: whe on 68 t 
H ‘ ° 
Me) a 
+ 
‘ 
: ‘ 
: . 
5 e 
‘ 
: ® 
: ; i 
: ° ‘ 
e ‘ 
: ‘ 
: ‘ 
: : 
Peas az 
; t 
‘ . 
a 
> Y) 
* 
* 
a 
. 
. 
: . 
e 
Ree meee cs cms See Sees FOan ee Otek eee oomasers t er ee ee re ee ere reer sy ~<— 
. 
‘ 
‘ 
a ‘ 
: a 
: . 
. 
e 
: ‘ 
: ‘ 
: ® 
° e 
: . 
ba oe - ee 
e 
: ‘ 
‘ 
° e 
: a 
2 ‘ 
: ‘ 
2 ‘ 
; * 
‘ 
6 
° 
. 
5 ‘ 
$ ‘ 
° ‘ 
b, ‘ 
=s - . ° sree} -~- ~ 
° 
. ‘ 
7 
a ‘ 
‘ i 7 i 
; ‘ 
: : : : 
. ° 
° 
‘ 
e : ° ° 
° 5 ° ‘ 
° 5 e 
: ; : 
* g ° * 
. * 4 ° 
Letter eee ee eee ee es P > $ : 
: seetecteos foe te Stet ts we trees 0 os ce twee cee me sccc aces 888 OF oS ORT OC OSs Co EEE SEEEE COSCO ES BBO HE SEs ¢ oe Sew b= oe TE Meee ee Be S* 8c eBees SAS SOteGere Prot r mee cccec ome toeeeces 
: ‘ : : 
e 
E ] ° 3 : 
: t . LJ 8 
° 
® ° i 
; ! i 
* 3 e $ 
' 3 . t 
: t 
H t : t 
! 
* 4 ! 
: ; : | 
We wt O82 BP © ote ae oe he. oes GR eo aece } 
eeort eeeeen - 
8 Ce ee © OO EN G+ o Sece = abe eo COST Chee Ob Staee as 4 PF octet ets wea tee ote eees cocemeceoettemscses OO we BE BOT Eo 8 OF OS SEH O Os OSE" SEH OS come co ef oc coe- erestors avs cases 
° 
$ Hi 
I | : 
e ® ‘ 
° Pq ‘ . « 
° 
; ; : i 
: ! ' : 
‘ 
a es 
. 8 = 
$ } : ; 
! « 
° 
s 
: | 


° 
Mpa nec ates pomecgroqpescccsom gros en sees inecs secon socnesenents euntenccaccenes o- oe Beno meee en come cans noe cc ccneee St cctens ecurece sss 


owe cece re ets esees 


Bet ee B= OBGe= .e8= 


19 11N0 BOLOY LY $3 TONY FONS IATL Vb a 


SI 





BI 


soueqstq TetrpeyY 


(seyout) 


98 





soTbuy MOTH eANnTosqy 3eT3INO 10370N 3A 
JO 3OTd BIO XTUORYaL pezyerousey oMuns W “ET sanbty 


(saeetbop) aThbuy MOTYA eANnTosqy 3eT3noO A0CIOY 


(UCD MOC cumt. “GOteiicCh OF Et (CoM Tis HES “6c %-Ge' 772 “OC. “G2 she 
OT 
a: : : : ‘ 





¢ > + 9 ~ “9 ~ is ; ° cl 
: : 
: : ; : 
om. ¢ +o ¢ - + misigsie ¢ oe 7s v ’ € | 
: : : a : 
: : : : : : 
; ; : : : 
rae — ‘ t +] -¢ eee t + eaceers “> . ? . p [ 
: : ; ; : : i 
: 4 e ry . e 
: ; A 5 H 
: : 4 : z 
t : : H : ; 
} : : . : 
: : : “ : 
genial haba deaba a oe eae se sae tnce neces —> “osces¢ os. onemecespseee G ] 
: : 
: 
: : ; : : : 
i i 
: ; : ; : : 3 
; : : : i : : 


eae fem pe mt eceeccpeeneecncntt pte fecstceeenee pose feseessenetdssee eed a pe tne ene Re ppmn ey aes SEDER Pct a! oe | 


. s ® ® 


a 

; : i 3 } : 

j ; i t : 

H ! t : : : : 

| , 

: : ; : 

nad cape coon! ’ sew atounes Tr Som ewer es { mee en ewes Soueaeces ease asans 7 secesccates t owe Como met + etteatecsce pce Sodasacvegccscesasae: ry $ 2 tac aaa oeabteas f J 
i : 

: i i : : : : 

i i 

pao nee aa ee Se ee ns so iceecseedeeeseeseeis arent J UNO YOL0N LV. S37ONY ADTs PuNTOSeY at 


eIJUeSAStTd (ervey 


(SauouUT) 





SOTbuy MOTI UOTRZeTASd 3OTANO AOVOY 3yA 
JO 3OTd BIO XTUORYSL pezyerouey omund VY “Pl sanbty 


8 eee ses Oe eee: Ly ace ee 


— 


° 
e 
’ 
° 
2 6 ams OF OF HE SERS OTe e 
. 


(soo0rbep) 


oo Os eee eee nets ame eene 


l 


er Pare 


6 
eee seeeef coe COO oO 0OF 00 Onn OP monde o Oe a ce O ee baam 8100. 
. 


Ome + CE ee Oe eE ret anetce 


22a Oo SOT se bastae 


2 Oe. C008 oe P88 8s Coase 


* 
99 meee th Cen eves et emes emer es es 


9 


@aaeee 


i 


>: woe Oe oe wooo ee ce 


8OR SOS CORFe Fs a teceaas 


2 8SP Ow eee sen aes Oe 8e 


cere pees on ces amesens coe 


aTbuy uotzetaeqd YeTANoO A0VoY 


S 7 3 4 i 0 





aD? 6 emt Shae seesetete -Broeess 


8 ee Obs ow or SHES EES? 


see cOete reese abate Sron 


= Cae HeFe ce erae SEne tease 


Porecceneccanes@eses seeeat. 


Soueeeeesersserneecees dt AMO B9108 37M HOR VIAIC-A— 1 9 1 


aoueqstqd [TetpeYy 


(seyout) 


100 


_ 





AQPOOUSA [PEXY Aeslul 7A99geAS Syd 71O 
4°OTd 8T9 XTUOFYSL_, payerouey OGUNL WV “ST eanbtg 


(99S/37) AQRTOOTSA TeTXY FeTUL 3703" 4S 
O6€ OBE Ofe o9€ OSE O7€ O€€ Oce OTE OCE 062 08c 






SOSe OBESE TEE I aoe 


sO e eet tt Ceeyecce Otero conert ta» 


. 
4 
SCO SS CES CO 6 OS SS OHS FS 680 FF 6 6 BEES © SO HOST © 

: : 

. ° 

* 

e 

ry 

* 

. 

+ 

J 

° 

a 

° 

° 

. 

J 

° 

1 

- 

‘ 

« 

es 

« 

¢ 

. 

e 

i] 

e 

° 

® 

CY 

e 

. 

~ 

. 

* 

es 

e 

° 

° 

° 

4 

e 

* 

e 

* 

rf 

° 

° 

. 

° 

. 

‘ 

° 

° 

° 

. 

iJ 

° 

ry 

* 

. 

e 

e 

. 

s 

. 

‘ 

. 

‘ 

° e 

. ° 

errr ee rt? tT eer Tire ee) e@eGeeeees ee :ete eO- SF er eebesGeoees 


ef we cc esse se eee OO! §- cece case beosbe-+o8 


eee Domne 0 cae 1800 wh Sows} euStarssce eteee 


as 

: 

i 

eee re TT tir Pre rr) 

! 

$ 

i 

t 

i 

i : 

ener eee cesses coc enpa eres 00s oo en eee ee 

d 

: 

i 

: 

eccaee ee ee 

i 

: 

+ 

ee 


© See eeegawenese 
ore 


tora de eeebeenesene ASNT BOWS ity ME sous! 21190934 iva Joy 


SoUPASTqC TeTpeY 


(soyouT) 


JO 





SeTbuy MOTTA eaAnTosqyY 3eTUI A0WReIS S3UIA 
JO 4O0Td 819 XTUORYSIL pszerAsusH OGUuNL W 


‘9T eanbtg 


(seezhep) aTbuy MOTYA sanTosqy 7eTUT 109e 3S 


Lt 9E ct Pe we a 2 ao 


: : 
: : : 
3 : : 
‘ . 
; : : 
: ‘ : : 
} : : : : $ 
3 : : ; i i 
s aia eels ty 3 t ® 
= aoa ont et aneee orn (at nial etananar Tenet: creer ceweerewen: seetce spt e@ec ssi nn teen wav wie «Pu cas wat wee soaoe: onan onh secccesanerenes matnere G te - emarteaceae 
: ' $ : H 3 H 
: ! ° : 3 f : 
‘ . H ° i 4 H 
: 
: : : ; : 
3 : : ; 
Wout ban oaen was oes Satan ge ges ose eas se aees bv gee ete ok cong as mee con ve Comer seam == les. -- ea@ee = 


® 
cf] 
® 
4 
i 
; 
eee 
| 
H 
i 





cece TPM BOLVAS 1Y $27 ONY ROT Ihosera—_ gy 


Souvistqd [etpeY 


(soyout) 


2 





SNTpey JFSETUI SA OTRSY SANSsoerag Te OJ, 103e3S 
284U32 JO 30Td B89 XTUORZHIL pezezsusey Ogun W ‘°L1 eanbta 


OTIeY SAnsserzd [e}IO], 7A07eqS 
GebG SUES creo Oo oo See ae OCS POCO Sell. TCSeES 8706" 

















3 
: , : awe 
: he ME eee : owe 
‘o- ? ® ' c BOG sce ® eo -¢ tcs6 : : 
° : ° : ° : : : ; ° : 
; ae ie ee: ere Fe 
"4 : os : sa sie bes PS eh: | 
: 3 one oe one oe an pe Gera src 4 Pes ae: ; 
° Mf . : : 4 : ‘ * H . ° H e : . . . C £ . 
ee : 3 io ; same deli) 215.5% ioakes ' co ee a 
6 : : ; i i a +f : ° . . : ‘ : : ° ° Ae a een Posey : . e : at 5 rs e a 
&- oor ‘: freee t Teer Pees SC eee nnn nee a sca a cpeet, Pepegorenct peek { ; : : : ' ! : : 
re Saker. Mee Meee phe So Race tris : : : eee es oe 
a ee es sa ye ane: eee : : : See 
: ° ° ‘ : . : ° ° : : : ‘ Z * : : : ei : ° : : : 
eet ae —_ ae _.- i a : a 
Sam eke Pars es bos ses tiene : ee : eels Wc tiueeceeas aa 
ae te i Wee ue ee gee a ves Seas ae ee ae : ee ey ee 
a : e ° e e ° . is Hy Ht ° . Py A 5 5 . e e * 
eS he ek ee eee ee eee eee oh se ee ; EA ae ee ne 
°, PS OUBEOLT HELL ROO E sor feces of -ceeepoe. sop oe Ce eee ee ee Pe | . | e t +e t- e-§- > FG de i a ee ~ ot t 3 te ee ‘ e 7 oo ' o. ! $ H 
: e = . . . e : * 3 ‘ . H . - rs . e r . . . r . ° : . . . * 
. Py + . e > ° ry 1 e s s Ls = - . =) : Ps r é ° : : 
: ee mec ea a artis wee SP Re ee ce Sone ee et a eee 2 Oe eee 
° ° . : . . s ° ° . : ‘ ° . . ® : ° : . e e . ¢ ° ® 0 
. . s . . ° . e . ® . 8 
: e ° ® t ® . : : : ry . . 0 e i. id L ry i. : . . 4 ® ® H : 
ee es > ee ee: eee tbe on Do ee ee he ies cae Pals ee 
° . . . ° ° . : A 7 Ps . Py ° ° ° . ¢ : s e : 
hs i ae er ee ps a a ae re ee erm ttt: 
Md bd s = e - e e e . e e : a 7 a Hy . e . . e 
: : ee : es , so : : : eS : : : ‘ : ‘ : 3 : : : f : ants 
A er ek Oe ee ee re eee, ee a: eas get ee eet eae 
° : : ° ° ° ° ; ; 2 : - ‘ : b aes eee recs aei@ce>.e emecae ay Ye oe te 
sonar ose: athe iicnesagons t a Pro-0:pocccep + Pets Pecce Ge o-Pee-- Peewee fecmonp =< e i mua aa 5 . : “f t + $ ! p ’ iv 
r) ‘ ° . 2 . . ‘4 . e e ® : t 5 H ° ° ry ° * ‘ i O R . : +4 . 0 
BOL fee Sigmon Ce Boke © Sets a ce ae ee ee Oe 1 es Se ee ee Bio ad 
ee ed Co ang i |e anes ee se es eee 6 ON A gc oe gee eee 
Se ann > bt tea gas aes : pe et ne eee to iy ier 3G pee ed ek Sh net or ee a eee 
So es eee es ae ; sage ieee eh ses ec een Venn PO Sa oe oe OR rc coer. ke 
: 4 3 i. H cS ° Py e ° e ° e ry Ps : ° ° id s = : 4 sl ry ° . e 8 . 8 
e e a ° ty e e . e . e e . e ‘ R . : . 5 ry 5 ry . e . . 2 . + ‘ 
: c : : : : : : : : ° : : : : ‘: : : ; : : Q : : : . : : : H : 
: - 5 4 : . ; . ° . t 5 ' e : ° ® , . ie ° ‘a . ‘ A . . ry rf ry . 
Secs ME sar otoe cnt : ee Ss te pert ey cote yer re eee De Te aan ne ite amen OR (MER ned amas 
rene Scene steer ee $n 8 - ge eGo ert ah-- oe >. omee Pome GH -ot pPoccea. Pree frets Focs sY ote etP Pee Becce wp + 0s me vane) Scar tte oe i aia 4 ' : ' t Ms $ : 
e ° e . e ° ° . . . . . . . . s Q ¢ : Le . $ x 5 . ry e ° ‘ ° . ¢ e 
H : H K ‘ ; jz : t : : : 5 ® e * ® . H : { ‘ ; . $ e H ® ‘ . 
H 4 : : : : H ° . A . . . ry : . i ° . 4 J t ; 8 ° ° ry rs e 
ee to ees ow te Ue Oe es ae ee ee ee ar eee 
: H H : ° ° ° | ° . ° e i . 4 . ° : a 2 ° e ¢ ° . : . 
: * : 5 : : : H : ° ° ° : : ° ‘ ‘ : ‘ : 4 F : : : ® ° ® 
® t e . r . : . . © Q 3 5 ° 4 e : ° ° . : ® e : . ’ s ¢ e 
° ® . e . ry ® e ® ' . : y 5 ig ry e . 
ei ce ae phe ee ae ae ee ee Dee Ca ee Wee Oe ster 
s . e e e e . es e e e . t e ed - e ° e ry e ry . s Li . e : 4 
CG aS Ge oe a ENE Ras eee ee ec cre SD eh ed ae Glebe ae Se cles lhe ec ek es 2 ey ie 
Save been fiteee geen bene cegeewee ges mahe er ert ceena wp raas eas ¥esiaiacde ae -ees gecceseu = eres or “gee 4 oor | “Prom p ree Rewes Greaes fences fescenfenceafeewret se t ' 7 ’ 
e ° . ° 6 . ° x ; : . , e . a e e z H ® e é e e . s 0 Q s 
eee ee ee ae eee ce Se nee aig ee a oe ee tee 
Dag pes SR nes cd 7 a see) ge a ee dae ae prt One ae ge ag ee ae ae ee 
. | . . 3 . . ° * ry ° e ° = rf t e . C . LA H rH H bd 3 3 t Py ° . ® 
. e e Py ry 4 e . e 4 . ry e . Py . ry ? C D : s : : e $ ry Py e . 
: ° : : t : : : ‘ ° H ° ° H i : : ; i : 3 ® : ° ° ° 8 . . $ H H H & t 
: * : . e . rf : e rf ' ° ® . ° . e 3 . $ e é . { . ‘ - a 5 q e t 8 . 8 
H . H t 5 1 . ‘ : : t : 5 4 ° ry . ° ‘ . s . 0 * : H Hy t . 
: . : ; e 3 i . H t e ® ‘ e 8 t ‘ e ° 
e e . e e i e . e . e . . H t { ‘ . i . ° ° e : § c R 3 
i : ¥ 4 : : H ‘ 3 : ‘ : : a ‘ e * . 8 ® : 4 $ 5 ® : ° ° ® 
eee ae a 8 a sic co Ge Cem eee OMe ame Maer tees AU PR gO ae a 
pererys Seat Soneny [neat Seneny Senne Dee Pros dsdee gem gees d ce gare: Rese fremeafercnepice tpetee eseepeccesptnes pete geet eb ect aets beseefectangitesefencengease eae ae ae 
Seer tia Gas veri pee ROOT eee eae a Ens Eile Sauna ag amma (A cane nage |e ae Se ee 
i oe i a a a oa ar ieee Pee (as coe: SE a ae Oe See re) ec cee ee ee 
+ . ry ® ° e ‘ e . e : e t ° ® t t o c $ . { . . : 4 e ° t i . : 8 . H 
; 3 ee et ee ee ee ee ee oe ee eee Ye ape cee eee terre gn eee emer ees A: 
ee ee Rete ee I ee a es ee oe ie ee re ae oot s) atatnint eGsa! Lis ollnieidlbc aneast 
. ® . ° * 4 r ° 8 ° 4 ° . 
° r a ° e Py ry ° rs rf * . ° * * ° ® ® t . : : i ° ry ® t ® 4 : 3 ! e r 
e e ® . ° ° ‘ : | . ® ® ry t : ® ¢ ¢ . ® 3 $ 8 2 x A 4 i i : : ‘ : 
. . ® ® . e ® ® e t : H H . e e . + iB . g 
OA ee ee oe Oe i eer hn en ee mee 
5 , ° . ° + e ° = Pf s ry ry : & Py : 
Oa ad iar a i i a sei aot tke ts wOpvis “Olive FwnsS3u0 W0j-ai 
- saiava saeimececiscanesee DS PP ee he er Pa en ee ne ee a Oa ee OR SN eI 00 Mates &- . Se ee ey e. CS i a Ces i ere oh 


SoURISTG TetTpeYy 
103 


(SououT) 





SPiPOCToO Ve [etsy  JoOTINOatoreIS Sua 
JO 34OTd 8T9 XTUORYSL pezersusy OMUNL WV “st sanhtg 


(9908/37) AALTSOOLTSA TeTXY FETANO A04eqQS 
(@3 


OSs O57 (000 “tes Cc cit Cs7 Cp Seo OLE OSE CSe OPE CEE Cee OTe Cor 66 - 
+11 
4. j ; 


—— 











: COeee-t otee Bie cieice ae : e f eoe ® ® Sn aiercien e a ' : : ~ A ad tee : cl 
3 : | 

qieccustbaea re cine Wei casita. ee Asko keer ere ‘ eke On ee. es e Be = - a i 

: : : : : . : fi . ; : H : ost : ete t= a . e. Sees s 3 ] 
: 

: : : : : ; : 

: : : 3 : : : 

: : : : : ; : ; : : : : : : : : : 

; i : 

: : : : : i : : 

: ae - ie - ST yor de en ae i e, ° ar bi oe< 68: eri - 9° seeisticlecaceesines iat ac a Gj . pcan “s elas vl 
: ; : 

: : 

; : : : 

: ; i : : 

: : : ; : H : : : ; : ; 

: : : : : : : : : : : : ! ° : 

: : : : : i { H : : 

° ® eens eccee- ie ipa tdi Raia sO. wm 0 00d wn 0+. cwee cet toca: a tatree me eee. sat PEER Peon Gaccsesiees Tore eee | oo q Ne eee eee ee e 

: - : : : 3 ; ‘ : : : : : : : ‘ St 
3 
et ee clon j.. ein ie Wer st eee? le | re 

; : ; : : : : : : ; 

i : : : 

; | ; 

| 

‘ nee wee sie cas otc msec st hescsessccnobeseneconetedjeesennsertafteoteesesse de eanccuceus | ee me passe oy eeteesteees DeCCCEEey St aPotccee o. hosts sceccteatoceses: a PACT harscs - 9 1 
: i ; : : : : : 

? ‘ : 3 H : . $ : H Py H ° 4 ; : 

7 4 3 : 4 t s 8 e a e e s H s 

: | : 

TO Geese: ame ge PAS circ, Nee 

: : : i i | { ; ! ; i H i : i 

+ e ® ; . a : ry Py j ry 5 H } 

| i } i : ; : : : : : : 

pies sii: a eee: 8) : i ; : 

: ° : - Bisteg ora Tipe ease daparcaiags Seer aa Ee aac cr ag euch Sud akiahicaes ooo ae a —ohny pemae shsceretaes l T 
i : : t i i i : : i ; 1 

| | | 

2 ee ee eee 

: | , & 

| | | : 

3 t +4 : i : i} e . s * e ° 

: et tlt ET, OTA Cet, rene Oe e2Seete on fcovsetecctecececct ee. | AIM pOLWES ty Fi soud fAir90TpA Wt xy. Ste 8 T 


SOUSASTdC Tetpey 
104 


(soyoutT) 





© OF eR Os 000s POF SHES OF OY OF 0 EES CY FON SHH HOBOS GSS BOTH EBLOSS SS OGOSS04TER- O20 


* 
{ 1 ‘ 
. 
pececccense antes et eee weacecesecccesd-ceccoce 
+ 
° 4 
* 
: f : 
fe cceccccccos rrr errrrrr? Tittre Tt itt. Ce 
e 
: : ° 
i . 
’ 
$e cce coqwwnnevetonsccs Proenevew- o-oo «G0 e ter Gone cnren ec ovececccceussecccoce 


+ 
° 
t i t : 
e ° . 
eer errr Perert fiir tte eoee Bsusdectenenes sas as< . oi . 


I : i 4 


ry 
Petree Set Poe Ti dere ece- 
+ 














Oe OO OTOL OOOT OS SOP ra Tr SEOF OS SF Es O* Or ges 
* 







a 
’ eeeeceerae eeoece rt @eeeat-eee @ « sa. @ 
Co mec ease--- 4s cena eraceeres as ' ° , - 2 « erete = ae 1 
2 J 
«eeee--crccce: : eee esecoocce Bdeecs soon: coerss « ond 
‘ t 








. 
e 

























OO ne ee ee ee 





























Peccceres cae comes 
. e ° | 
: ; : wT 
° ' . 
Wererreeerrirretrorrs oor Steere Bp cnes- te ororeeyes « Oc emereeeseeseces ec esenseeses a eee Pucccovemenses: ercee A 
rf . ° 
{ ie ry e . 
* . e 
: s : —t 
oO wwe Ooo OSs FTRO SS OOSS POSOSSeWS CBSE EEET OF a@eeee *20+ O08 ser freceree eer @ +e ° oe fe Oae Fes che sarsegr erases sseeeseteses-« 1 
. . 
: ° ‘ . 
* eo e . 
. ld ° e 
oO Owe cow SOs Serer ousee PO 0 oe 40 00454204 Oe EES BD 0+ O° 5 SHTES HOH OFOST GETS S- oe ses ae - erro Sei OI POOC ert t er sew 
e . 
i MH Q . e 
Py I e e ° 
Ow tee SH es OF See BH eHOS: BS OHHESSOCOR S06 60 ans Gre es's2eer00" See cece a* Rpt TSK SETH 2K Ee) . SP OCC eee OF FO 88 OSs FOSS OS FOES SEs se HR000 
. = e e 5 
: : : ‘ e 
: : e : : fe @] 
Prose cwwwcemees emenese-p- cones eee ee Ter ee er eo e 
* ° 
° e . y 
: ; i : ot 
Do - 00 oc eae ew ceemeses caupacess Groce 00 wt ence: Goes aseeee seawcs ve cees de cces + cess cnesecs crass forte ccacowesesesecees fos-0 ee 
e a 2 t 
e o ° ° 
i e s a 
esuiisscercetcenteme esp sce Peeteteors eeneee PUT eer Perr © Cec sogeocescecomcscesse.ces | cHe-sae cones amen “As 
e ri . ‘ "5 
. s t « H : 
° ! : : : : ' 
b. eee-ccee emeet0 wees co ae 0 Ge 80 Ca eha ec O0 OF Oe 2Oe Es He FES) C8OSS- HOCSH TO OVO 8 e © 8 OF Soe CeO ORes BOROEs SeFeBEnEes Se «: Mth the Se 5 
e * . * 
4 { e 
a e 
° e + = 
OS Ll le etd J CODS COCemD OB EOS COS + H O59 00-2 0008 C8002 OL ODOM Ee oe . cos eocetoocesces cone seee.re "eee Silas ses Sr 
e . t e 
* e e 
. e - e e : 
oe Oe OSe OH Os Fe HOE CO" Oe Eh OOO IOS | eos sere. es 08 Gece. eeterse-0e © fe Jecwcss cccwerescceces +o 
e e ° 
! : : 
* * 
ene. oo See tuaneceos 06 ume ewonceccc cee coe Gaweeraesee-se*es & -4 ere o-¢- eo Crweses secme Serer esecacesese esse 
e * ° 
t e 





e . 














« 
e 
e 
° 
° 
e 
° 
i 
. ° . ° 
$0 we. tees ecesconwase oc cee cut rene secon cress quenccecccen ccosecces: é- CeO 8 8 OH ee BBO! CEOs OO OS Hes Bsc errs HOt eees OFess eres 
e ‘ ’ * 
: : . : a 
° . . W : e 
Ge -caecccaracereree Peretti ti rt eee i er oer eee ee Ter TT eT Te tT Perey Terre rr eT? tt) tt) eee ee ee 
e e e C e ‘ 
e e i t e 
. J e ’ 
° e ° . a 
Pc cece coneeve coeme © 06} ott ow- ames coeessoes o- Geo-s =e. se te ° ¢ t * ° eF.e8 es80 ° o-. e- eee 
E . 
ry ry ® ry ° . 
4 . 
$ H : . ° : cS 
P wwe oe ercwe moncee 2° © woh c Came 00 00 HOR CaOT Oe Ss Hew eres. -eessacsovor-« ¢ ee Corre § Oe er - cence re erereseP: oc: aerate se eacersess s e 
e t + . e e 
e a LJ e 
: ‘Ss e ao ° ’ Na 
[ et oem ah +6 8c cece es comcame ees cu ee teers: © eee eee er errr eres Cre ere es rere o§ wwe swe eee cocesracesss ' 
e . . . 
+ i e * e ‘ 
3 ry ° e e + 
[} Popo IS whew Cone ce sconces cease ceees: Gesacce-se ce secese «ey * wermerce +s fececcecaeeesasce © eeecs ees ees eelenlcciciasis oe Cc” 
- q 
: i : : 3 
e e e e e 1 
7 e bd bd e . 
BACH ol Oe 6 Pe Oc Cle £6 O0E TOO SOO TOSS $ TOST OMS OF amEO+ OF BEDFe+s eeEL ance ee eee ee eer Peer rey err ei rerer Sree ry) 
° « e 
° ‘ | ’ e 
: ! ° : : 
° ’ « . . - . 
Joc ecceccee: Pee ee ee eee ee Boers eco ee ce ce cee mmecseces wee Coes ees ee ees Oc Os cae eer-Oeee e-* ee 
. * 
. « 4 e e 
° e e . 
e e e e . e 
pecs ceccece POSS SSEES SP O00 FERS CFC EES C2 CODE SMM Cee- S555 0-0 "es = "eevee 8802 @Orcee t © eee ernees Or © eaf cer esecccetearecs conertooe eeere 


= re sow sense Haceeecomneerenmmenegioees ores worse =e te vse e+ 






ee ee 


- 8H Cw oees Coeee he coees Steer aeees- oes - homer cseeeses-scascaseshoaessn sees 


‘ °° 
‘ ‘ . ° : 
ad A ° 5 e 
’ Gavece-8 (enéaeaeoe Owe 00s 00m 000 Ce ew ee Mh aameewereses-cecse 2 ay bl oe “8 finwes sibs ees vcoeeneevemecceunvesen « 6 ce . ; ee ! 
‘ ° ‘ : 
2 ‘ : 
5 ° ; 4 
5 : © 5 
@ seerse © ORO SSO! OO CWO! OS 09 08 OF HOD CF 0 GREG OHO WHET OrEs -B O08. 4 eee eee er eee eed . 4 A 
. i} . ° e 
e e e . 
e e e 
° : 
aigach aaah scasecasanetcousanecedbeotsclene:beoueas6 er. e eeeeer Sarees et *eneOe «ee 1 e- e. 








° . * 
: : : c< 
eorce Ga. 20a mors Pesemes ccceceewoenseve Pres sceeee-.+s = ©@e + @ + eve Pe Oeh ees FOr Ors car Ose- 60s OO Ose Pees ae -seteretesecsese 9 ° e 
5 * = 
fon : : : oe) 
teres Oe eee 0209 000 ap oom en es oF Hewes Sete: cow Pe ecoeccoess-s- -s see hwe- 0-2 --20 20 eree-§en o-0 + or- oo Cre Ce eee ees Pee ' 
: : : : 
c q ° ° 
SW ce PS FOODS © O% Bint D+ FO ORO ©0066 OF OOF OED O44 HTH OO O4FD 5 Ot SOFH EH ©" - errTirrrirt Tr ttt tT te rr or or Sr See OO - Aree see Orr ne anrs nee “Jj 
od : : : 
‘ 8 H = 1 
7 see Gee POSS OSs 0 O6 6D 66S 6 ONT SS OS OOS PERS OS OE! C8 OF OO FO OF Gs FORE Hs O* OH EGE: = OOO OOe oe Oe: Cee 006 oe EROS Oe OOOO Foc e Ores BESO FO eet es +s Os Or eee =e @f* a2- 88-6 wore 
‘z ; : : 
.< : : : 
O + eeer-gae OT Oe Om O00 Oh Oe © OOF Oe Oe OOF SENOS SE CSOs OOF 200 HOES OSSO Ss FE- 8 UE EES ST EFS SHEE OL FEY HOSE OF SORE SS OPES: OF SEO Ors smBEEEDEESEI Es Os 6% © ewes eee ene 
e . e 
e e * 





‘ . 

Ce ee ee 
. . ’ 
r e 

' : ; 


ee et 2 Poe Peer ees er 
Fi 5 = 


FE FUG 


eiantanauicsssecpassacsen sees eeasieen se qussescoseses everserscas 


t 
> : : 
4 






> Of cc wc OOP Sewers Ws a OSs 06 02 OSES OMe COTEEE OME HOSSS SS areerers 28 forenee 





18 
17 


oO uy 
— 


(Seyout) 


° e 
Pi da ry ° q 
‘ « * 
6 Co eS 6 OSTT TS Ee F OS OF 00 OF OF 6 CSREES ESET © OF 00 O64 0088 0800 08005 86 88 OOF 00 OF EE OOSe OOOO HOSS HOOF OO TEs OOO OS Os OOOO wee 40 Re es 82+ Hak cee mh erereeeereees see t 
. = e . e 
ry . 
. e e 
. ’ ° . . 
megl> «§ coc co camemmszeece POO Oe 0 OOF 2 OOHRS OS FEF SPAT TSS OF OWS SODSS- Os0% Oem - Oc eat-+ seers steessese Sots ceete ese Coe esesemesers- 800 «- . on fet ee + Omer eene a seee 
w- J . 
’ e ‘ 
= 


SFO OF 6 COO FO OEE FOO O08} CHOOSES OOS O ESE Ee FEE A OO EOe Serer eres esas 


"eere-etec#ero- Pocwe covccocscceccccene§ cece cscecccccecccasss ’ 







ee ee ee rey Ce ee ee ee ee eee 










lawnelwelcesisienGecns'novislaeicein¢) Suiew en) Plesieviesulsaciisiesieseiesis) § Cesescwe x as. 4 oe 






a 
a 


13 
4 


Sodue4Stq [Tetpey 


uO> 


(degrees) 


Stator Outlet Absolute Flow Angle 


618 Plot of 


A TURBO Generated Tektonix 


oe 


Pvqure 


the Stator Outlet Absolute Flow Angles 





SeTbuy MOTY vuoTQ_etTASed 4eTIANO 103e3¢g 
94} JO 3OTd BIO XTUORYaL payerauey OGUNL W 


(s90zbap) aTbuy uotyzetaAsg 3eT3noO 103e 3S 
O{ O76. 956 0) 2°65 5 g°e 9° pee 2'R g 





ec: : ae: : Se 
cane: a: na! 
: : : ‘ : : : > ; i 
: : : : : : : : : ; : 
Bowes Ge. = oo a oe a. : 
nee oe Sac ae es Sr ts: ee aes : 
es ee ses nc cernon i Set kt es rs 
ee eee ae ee ee ee ee ee ee 
(Se a es ee een 
ee ah Oe, ea: ae. a oe ne : 
ae eee Sec mp = ot s ° H APSO -~ atc e cms en2e aacee o 2 & owe » oe AG me abetecee Scere Cm “te nies. 
Sab alate oe” One. che Ue ee en a ee er eee 
i ee fe Gg Ge ; 2 ee ee ae Be 
ee el et ee er, 
ee oe ee ee: a ee ae ee ee ee ee 
i pees l i omaut : eo etd eet ie eee 
. . . e . $ : +4 . ; { : . H 
: ; boas ie ie eee ee 6 SRR ae ie eee ele 
eet ooenee Caner ea 4 mead ieeees tee oes tee Seaton Soe Deane a vt: G-ctmrebeceneede cone aoe ie peers cs i ~-4- See e 
eee ah ns ct mu fos hte Ati Te 
St gece gy cs inst tages a, nc aes SSNS ee op: * i 
oe cee ee ee ee ee se ee ee oe 
fo Re eee es: ee ee ee ee. Oe re aye 
: . e . ° ® 6 * e ® ® . ° ° : A : : ° ° 

° e ry ry e ry ° € ry e ry e rt Py 8 4 e $ ° i 
OY "2S NERS et oo 9H) ita, Sia ne eg aC Se SS ec 
een ee eee | a oo ee ee ee 
-aadeiet Semmes Senne entans Sinnaah caeatad Sani Siem cates Sieetth Sotcels eteeias Santis Coneaas Cote Sete eee Cantor Coton earns Unies neo 
‘OY CC yh ree ee et oi ~ fe Coe ste yg ae 
a9: ot Ciera: we a ee es 

s ' ° 4 r ° 

a ee eee ee eee ee es ee ee a ae 
Pie tt Se es a ee eh a ee ee 
iat 2: Tt Pb 
° ¢ e { ¢ { i ° t H ° . . t . H Py ° 
i See tat : | : | : Sere : : | : 

4 $ e Hy H t e ° . ° ° ° ° 

Hae, ¢ ; bangs! eet | : ; 
eae ae os Ae eae ame ea ee ’ Seal teS eats seat es Cae Le oe ao en begs Soret 
inet i ae | ne: 
a eee : ee ae a et ne 
ee ae i Tipe a ee hes A | ee tte et ae a 
| $ { t 8 | ° } { 3 ; * : ‘ : 8 : ° 8 i 
ee ee es Se (ee metas Tn Ee re ae ae a 
eee ne en ee et a a a a oe ae 
i oa es : SRN eet PC il t i ; 
: ¢ t ' t t t ° $ e : : 5 5 ‘ t 

“4 s $ . . e $ . ° } . ° s 
en, eee ESER Snot eee Ae eae, EON eee B etoeae ie Some ® a a Oe an Orns oo pane een SOs eese oc eee AY 


Someees tate peehsteans. ce -a-e 


"OC. 2 nb 


Bee? ee, cD 


im 









EGR ale a 4c | 
Mercere 
ee 
oS Se 
ope ae 
iat ie 
oa 
ebee beg P. bp] 
a 
he ae 
i ft 





| 


e 
s 





i 
! 
eernee oe 220 fro coe one 


SSebe oO 20 CHESS cS 


4 


ASM: voLf1s PMY HOLA — | oy 


(SSyuoUuT) soUuUeySTq TetTpeY 


ING 





ubtsaq 30g ‘ARTOOTSA TRTXY FETUI AOAOY AT 
SUOTJPAATISGO puke OFTASHO OF SUOTROTpeIg Jo uostaedwopj “Tz sAanbty 


ALIDOVWA IWIXY 


“OSh “OOn “OSE “OCE “0S2 “002 ‘ost 





‘9 
“Ol 
“21 
za 
nD 
oO 
ie) 
ba) 
| Sen 
‘at oO 
wn 
—_ 
wD 
ra) 
rm 
‘91 
“9! 
OLJAMO 16 S21NSIW r 
SIN WA 03131030 F 
SINIWA OFANISAD @ 
9°40! = *O1y SSM 
N91690 20S © 0334S 
"02 


LAS INI YOLOY 


107 





upTSog 205 “AJLOOT2A TeITxXy 497139n0 10704 107 
SUOTRPATOSAQ pure OFTALH OF SUOTROTpSeAG FO uostaedwop *77Z eanbty 


LLTIOVWA WIXY 





“009 “00S ‘00h “OOE “O02 
‘Q 
“Ol 
“21 
DD 
»D 
oOo 
me) 
on 
‘ht O 
| 
72) 
“+ 
2D 
z 
Mm 
m 
“OT 
"el 
OJLAWS 19 S1INS3U 
SAN WA 0313103Ud Y 
SINWA O3AN3SE0 © 
9°201 = HOVs SSW 
M91S30 20S = 0334S 
‘02 


L4JVLNO YOLOY 


108 





ubtsaeq 30S ‘ARTOOTSA TRTXY RYSTUI 3A0ORFeAS AOJ 
SUOTJPATOSGO puke ORTARH OF SUOTIOTpetg Jo uostazaedwog ‘Ez sanhty 


ALIBOVWA WWIxXd 





"S2h "COh “Ste ‘OSE rd 3 “O0€ *6L2 ‘0S2 r 
“Ol 
2} 
D 
»D 
oO 
D 
= 
“ht oO 
w 
—4 
D 
= 
ran) 
mM 
91 
‘es 
OL1javo 16 SiWS3u Mm 
SINWA 03191030d ¥ 
S3INIWA 03003590 OD 
e°20t © MO14 SSuM 
W91930 Z0S = 03346 
‘02 


IPE INP SONGS 


109 





ubtseq 30G ‘AITOOTAA TeTXY F#ETINO 3209"9S A0J 
SUOTJPATZSGO pue OFTARD OF SUOTROTpeAg JO uostzeduopjn “Pz aanbry 


ALETIOIWZA WIxd 


“OOn “OSE “O9€ ‘One ‘026 “00s 092 





‘Ol 
"21 
»D 
D 
oO 
f=) 
xD 
— 
‘al 6 
wi 
=< 
D 
z 
i | 
™ 
91 
“et 
OLNTAHS 148 S17NS3uU & 
SINWA 03199030 
S3INWA OIAUISAO @ 
9°401 = MON4 SSHM 
91630 YOS = 0334S 
‘02 


LITNO YOLVIS 


pad 





ubtseq 0G ‘’AQRTOOTSA TeTXY JeTUL AOROY AOJ 
SUOT}JEATESGO puke YOSATH OF SUOTRIOTpeAd JO uostzedwuojD *GZ oinbt4q 


(93S/13) ALIIOVZA WWIXB 


“Osh “OOR “Ost “O0€ "0S2 “002 “ost 





Ol 
mt) 
D 
21 = 
pe) 
on 
Oo 
wn 
~ 
i D 
nl z 
7) 
Mm 
ae 
m 
2 & 
‘or om 
hel 
‘el 
HISUIM 2@ SIINGIV M 
CANWA OF13103Us VY 
S3NWA 0FAN39SA0 O 
8°20! © HOTS SSuM 
#91830 20S = 0339d6 
“02 


LAIINI YOLOY 


aa 





ubtsoqd 30G ‘ARTOOTSA TeTXY FETFNO AORVOY AOJ 
SUOTJEATOSGO puke YOSATH OF SUOTROTpeTg Jo uostzAedwopD *97 SANbTY 


(33S/13) LLIJOIWA WIX¥ 


*009 "00S “POR “Oot “002 


“O! 


2) 


“ht 


JINYULSIO IHIObY 


(SSHINT) 





HISVIN 19 C1IINGIV M 
SINWA O0319103Ud & 
SINIWA ORANRSAD DO 


O°201 = MOTs CSUN 
W9IE30 20S = 0393906 


“02 


Hee WOMB slay ible! 


Es 





ubtsoq %9G ‘SaTbuy MOTI SATRETSY FOTINO AORFOY AOZ 
SUOT}JEATOSAO pue YOSATH OF SUOTRIOTpsaIG Jo uostzedwop °*f7 eAaNnhbty 


S2°09Nb MOTs SJALitb ted 


‘oc! “08 “09 “On “02 0 





0! 
D 
D 
‘21 OO 
D 
| am 
QO 
Ww 
ad 
4 TD 
nt Zz 
™M 
mr 
= 
(ae) 
a 
1 4 
Ay 
“Ot 
HISUIN 16 SI1NSIU M 
SINWA O31,/03e & 
SIN WA OFANISED O 
O°201 © MOTs SSUM 
N91S390 20S * 03346 
‘02 


13 tine YyOIOd 


i 





ubtsoqd 30¢ ‘SaTbuy MOTA e4ynToOsSqy 32eT3NO AOAC AOSJ 
SUOTRPATISAO puke YOSATH OF SUOTROTperG Jo uostaedwoyjn *gzZ saANbtyY 


S31VING MOVs 34N 10S8Y 


‘00! "00 03 "On “02 0 


“O} 


“eN 


Ted fOb” 


“ht 


ZINYISIO 


“gt 


{SSHINT) 





"et 
HISUIM 19 S2INSSY ® 

SINWA O319103d 

S3NWA O3AH2SHD © 

9°401 © MOTs SSH 

W9ISIO LOS = 03905 


02 


13 ILNO YOLOY 


114 





ubtsoqg 30S ‘AQRTOOTSA [TRTXW FOTUT r0OReIAS AOJ 
SUOTJeATESGO pue YOSATH OF SUOTIOTparg JO uostazeduoyn °6z7 sanhty 


(93S/13} ALIIJOIZA WIXs 





"Sen “OOK “sce “OSE *$ze€ “one "St2 “062 
9 
01 
» 
D 
in 
D 
[ io 
oO 
Ww 
a4 
, D 
a 
™ 
m 
z= 
™ 
¥ 
91 OM 
wall 
‘el 
HISVIH 16 CIINS3IW R 
SINWA O313103ud 
SINWaA O3AH3SA0 © 
9° 401 © MOTs SCUM 
N9IE9O XOS = 033906 
"02 


LA INI YOLEIS 


le, 





ubtsoq 30G ‘AQRTOOTSA TeTXY 2eETANO 10ReAS AOZJ 
SUOTRJRAITOSGO pue YOSATH OF SUOTAOTperg Jo uostazaeduoj ‘gE Bsanhty 





\ 
(33S/14) ALIDOIZA IWIXY 
"OZh “OOn “OE “O9E “Ohne “o2¢ “OoEg “0e2 _ 
"01 
xD 
mo) 
"21 iS 
»D 
i 
Oo 
wn 
—4 
: DD 
nr > 
-Q 
Mm 
z 
[ ) 
3 
“oO! m 
ek 
“O1 
HISVIN 14 SiwmS3y mM 
SINIWA OF13I0Id FY 
SIN WA OB3AUISHO @ 
9°20! = wO1d CSUR 
N9IS30 Z0S = 03346 
‘02 


13140 BOIS 


116 





ubtsoaq 


"OS" 


ie) 


OS 


‘VZOTUI TORZOY ou. Fe AQTOOTSA TeTXV 


“OOh 


SINWA O3151034d VY 
SINWA O3A43SE0 @ 
2°ntl = MOVs SSUN 

w91S390 20S © 0334S 


(235/135), LEV IO13A IITA 


“OSE “OOE 


LI INI YOLOY 


"082 


"002 





“TC. ounbD ia 


“OS! 


‘Ol 


‘21 


"hl 


“91 


‘ol 


"02 


“B10 


Neisio 


a2 


(SSHINT} 


Jae 





ubtseq 0G ‘3OETUI AORZOY We SaTHhuy MOTA SATIAeTSY “°zE eanbty 


SSIONY M013 SATIY13Y 





‘001 ‘00 ‘09 th.) “02 0 . 
“Ol 
mt) 
ma) 
‘21 OQ 
»D 
r 
oO 
wn 
— 
. D 
ht = 
a) 
m 
z 
om 
T 
rm 
91 Wn 
“ol 
S3INTHA 03191030d Vv 
SINWA O3AN3ISE0 O 
é°htit © MOTs SSUM 
NSOIEIO LOS © 03306 
*02 


LIINI YOLOY 


118 





ubtseqd Z0G ‘AORVOY SYR JO OTReY SaAnsseag [Te}OL “ECE VaANhHry 


O}1bY JYNSS3SUd WIOL 


s*1 fh'l i | ai | | O°} 





9 
‘01 
D 
D 
‘21 O 
ta 
D 
(= 
OQ 
wn 
—d 
nie 
m 
m 
Ya 
om 
ek: 
‘or om 
hs 
"Q) 
SINTUA 0319103Ud + 
SINWA O3AN3IS80 © 
Chill © MOTs SSHW 
N91S30 YOS = 0394S 
“02 


IS tN! gdbadoe 


1S, 





ubtTSeq 0G ‘1030N eUug JO ADUeTOTIJA OTZeqeTpY “PE eANHTY 


YOLOU ‘LINIIIIAIZ JLIYVEHICY 


ot 6°0 8-0 t'0 9°60 








9 
0) 
SINWA 0319103Ud ¥ 
SIN IA OIAUISAO O 
cont) = MOTs SSHM = 
NOISIO XOS = 0324S ‘24.0 
4 
D 
= 
s) 
wn 
4 
ee 
Z 
m™m 
z 
m 
2 
‘gt OM 
wi 
7) 
“02 


LA INIT YOL0Y 


2.0 





ubtseqd $0G ‘3eET3NO AORZOY syW Fe AATOOTSA Tetxw 


SIN WA O319103Ud 
S3IVWA O3AN3Sa0 © 
2°nit © MOTs She 

N9IS30 20S = 0334S 


(DASZLI) AILIIOWA WIXY 


“OOF "OSE 


hapa re stay ft gts 


“OOF 


"OS2 


-ecev oanbia 


"002 





01 


“21 


“hh 


“91 


02 


JINBLiS(IO IdslOeY 


(S3HINT) 


iene 


-. 





ubtseq 0G ‘3EeTINO AO-OY 3e SoOpThbuw MOTA sATReTeyY “*“9E BANHtTY 


S3JIONY NOS JATLN13Y 





"001 "08 “09 “On ‘02 0 
. am: 
Ot 
» 
1) 
| Oo 
D 
| 
oO 
tn 
4 
: D 
hr oS 
() 
mm 
> 
q , 
» 
gl mH 
al 
‘@1 
CINWA OF19103Ud Vv 
SINHA OFANISBO O 
L°W11 © MOTs SSUW 
W91S30 20S = 03306 
‘02 


Ieee ss 


2 2 





ubtsoq 0G ‘J3eETINO AORZOY 3e soeThuy MOTA saANTOSquy ‘LE sANHTY 


SJTINY MOTS J1iNIOSSy 


‘001 “08 “09 “on °02 0 





“O1 
D 
D 
‘21.0 
D 
Lee 
s 
wm 
oe | 
at oS 
z 
ty 
z 
z 
in 
Br te 
91 
SINWA 0319103U8 
S3NWA O3AUISGO O 
ERTL © OTs SUM 
N9ISIO XOS = 093486 
‘02 


f4alind Yoru 


123 





ubtseq 0G ‘3eTINO A0R0Y Ae SOTHuy MOTYA uoTRIeTASeq “gE saNnbtTy 


S3TONY NOTIYIA3Z0 


“Se Ce "$3 “oO! aS ‘0 





"oO! 
D 
D 
"21 [ea] 
D 
- 
Oo 
wr 
=~ 
ot 
™ 
Lae) 
z 
Oo 
ale 
"9! mM 
Th) 
*e} 
SINWA 0313103Nd Y 
S3INWA OFAHISAO @ 
dZ’ntt = wold SSUW 
N9I1S3Q ZOS = 033d¢$ 
‘o2 


L3tiNO YO1OY 


124 


: 








ubtsed $0G “IO8TUL tOVeIAS |BYQ Ae ARTOOTSA TeTxy “6E sanbtTg 


(ose724). SPC Ths Wel xXe 


“OOS “ocs "oc "ose "ONE 0S? O02 
|S aga a ae DE ena g 
or 
mee e : 
e 
my 
r 
ral) ena 
D 
| oe 
ay 
ur 
m= 
o ram 
L. 
ae 
lee : 
} “91 a 
a ur 
r) 
* 
“@) 


SINWA O3L3IC3Ibd 
S30 Wa O3AH3590 © 


c°HIT = MO1I SStM 
WSIS30 20S = 9334S 


"02 


Ire INO ete Benois 


125 





ubtsed 0G ‘3eTUI 2z09e9S Je saTbuy MOTY syNToOsSqY °*OF Sinbt 4 


S319NY MOTS JLNIOSEY 


“00! "ee ‘09 "Ch ‘02 ‘0 





‘Ol 
pe a) 
Dd 
“et = 
=x 
ts 
| ie 
w 
= 
5 D 
fl Zz 
> 
rt 
Pre 
™ 
Ag 
‘OI m 
in 
‘@!l 
SIN THA 03131030d 7 
SINWA OIAUICAD © ; 
Z2°HLT & MOT4d SSUW 
N9IS30 X06 © 033988 
‘02 


LAIINI YOLULS 


126 





ubtseqd 0G “‘FEeTUI A0ORVeAS We OTIeEY oANSSeTg TeIOL ‘Tp eanbtg 


OTLWY SYUNSS3Hd WI1Ol 


ol so‘t 


Td TOBY 


NB1Si0 


je 


(Sse) 


S3NWA 0313103Ud & 
= SINTWA OFANISHO @ 
2°nit = NOTTS SSUM 
W9IS30 20S © 03306 





ILS INI YOLbI1S 


ie] 





uBbtsed 30S ‘3eT3NO AOReAS |YyA ye AQTOOTAA TeTXY ‘Cp eANhHTY 


(J3S/145) ALIIOIZA WX 





‘005 ‘OSh “OO "OSE “00€ "082 "002 
x) 
"ol 
a) 
D 
‘et © 
»D 
~~ 
o 
uw 
ome | 
; D 
hoo 
m4 
ima) 
z= 
™ 
SE 
opm 
Wy 
"ot 
SINHA 0F19103Hd ¥ 
S3N WA OZAUISeD O 
Lh) = MOTs SSHM 
N9IS3O KOS © 0334S 
‘02 


eee SC ee pletl sles, 


EAS 





ubtsed 305 ‘39T3ANO 103e9S 3e soTbuy MOTA SANTOSqNY ‘EF VAINbTY 


S3JI9NU MOTZ 31N10SE8u 


"06 ‘02 “at “0 “at- “Oe- “O6- 





‘Or 
2D 
D 
“2h S 
D 
[ game 
oO 
Ww 
— | 
: D 
a 
™ 
my 
= 
Mm 
mi 
1 V 
*@1 
CINIWA O319103Ud F 
S3INIWA O3AH3ISeO © 
2°hit = mOl4s SSYM 
N91S30 20S = 03346 
"O02? 


ISI AUG) eshalial ES 


a2 





ubtTSeq 0S ‘38T3NO 3707e9S Ae SaTbuy MOTY UOTAeTADG “PF aainbht 47 


SJTINY NOTIYTAIR0 





‘$2 “02 “St “oO! °S ‘0 
’ “qQ 
‘oO! 
yD 
2D 
‘21 0 
D 
tf 
| 2) 
tn 
4 
. Db 
wl = 
(>) 
mH 
a 
i) 
aye 
>| ae 
ad 
"Ql 
CINWA O3913103Ud FY 
SIN WA OFJANIS8O O 
2°hUl | MOVs SSuM 
NOIS3O 20S = 033d¢ 
“02 


fe asta alnicialeS 


ie 0 





ubtseqd 3024 


“009 


‘ADSTUL AOVOY suq ye ARTOOTSA [TeETXY 


“Oss 


(J9S/15} ALIIJOIZA WWIXb 


“00S “OS 


“Oon 


“Oss 


Cr Sam Oe 


‘oot 





SINWA 0319103ud F 
SINIWA OFANISEO @ 
SS'1S1 = mOty SSHW 
N91S390 x04 = 033dS 


f4 INI YOLOY 


Ol 


‘el 


‘hl 


“gt 


“oO! 


“02 


INBLSTQ IWWidd" 


= 
3 


(SSHINT) 


Se ail 








ubtseq 3O0L ‘3eTUI A0OFOY Qe SOTbuy MOTA SATIeTSY °OF oInbty 


SJTONG MOTY JATLYTSY 





‘001 "08 “09 “On “O2 Le 
‘Ol 
BH) 
D 
‘21.«0 
D 
o 
wn 
~ 
es 
a) 
mM 
~ 
2 
™ 
as 
‘arom 
wn 
‘ot 
SINWA 0319102Ud 
SINTYA OIAUISEO @ 
SS'IST = MOTs SsuM 
NDISIO YOL = 0334dS 
"02 


13 INI YOLOY 


ro2 





ubtsod 304 





‘10}7OY 3sUuR JO OTRPY SAnsserzgd [e}oL 


SINIWA 0313103Hd VY 
SINIWA O3AH3ISAO @ 
SS‘°1S1 = MOTs SEHM 
N9IS390 ZOZ = 0394S 


OFIYY JYNSSIWd WIOl 


€°t 2"1 


JA INI YOY 


“Lp oanbty 


“ot 


21 


“ht 


“oi 


‘02 


JINELSIO TWICsY 


(SSHINT) 


ies 





ubtseq 30, ‘arz0R70y aYyQ JO AOUSTOTIFY OTReEqeTpyY ‘Sp saAanhtTy 


YOLOY “LINZIIIIIS3 JLLVBYIOW 


o't 8°0 8°0 £°0 9°0 





SINIWA OFLIIOIUd FY 
CINWA O3AH3IS80 @ 
SS‘°1S1 = MO14 SSHH 
NOIS3IO xyOd © 0934S 


13 Ne eae 


“oT 


"2 


“hI 


“02 


JINXLSTIO Wrasse 


(SSHINT) 


134 





ubtsed 30L ‘3ETINO AOFON SYR Fe APTOOTOA TeTXY 


“009 


“OSS 


SINTWA OILIIOIWd Y 
SIN WA OFAH3SAO @ 
SS°1IST © WO1d SSUW 
W9IS390 LOZ = 0334S 


(J3S7 1S PETIO VA Wixe 


"00S “Osh 


L37LNG YOLOY 


“00h 


“Ose 





“6% orNbta 


“oce 


‘Ol 


4 


“nt 


"et 


°02 


JINSLSIQ WIidkY 


(SSHINT) 


13> 





ubtTseqg 30L 432e8TANO AOQOY 3e SEeTHuy MOTA SATReTSOMY “OG sANnhbty 


S3TONb MOTI SAELY13U 





*oO! "ae "09 “On “oe 0 - 
“al 
pe 4] 
xD 
"2! = 
»D 
r- 
(a) 
—) 
Ww 
-~4 
p D 
Ln > 
o> 
m 
ra 
Cy 
a 
“91 my 
2 
‘et 
S3INIWA 031910FUd F 
SINHA OFAUISeO O 
SS°tsSt © MOVs SSHW 
N9I6390 202 = 033468 
‘Oe 


LS TLNG YOLOY 


136 





ubtsed 304 ‘#8TINO AOJOU 4e seTbuy MOTA eaNnTOSqy “1G eanbtd 


SSIINe MOIS S10105 6G 


‘O01 “08 °09 “On Oe ‘0 5 
| PE ES SD 6) ED ET, (ES ea aa = 





‘Or 
w 
D 
‘2t OQ 
DoD 
tr 
QO 
w 
=| 
: D 
fl Pa 
(4 
mr 
Pd 
Oo 
2 a 
‘groom 
wn 
‘Q1 
GCINWA O313103Ud F 
SINWA OFANIS9O O 
SS°1S!} = MOIS SSUH 
N51S30 204 * 0330S 
‘Ne 


L31iN@ YWO10Y 





ubtseq 304 ‘32eTINO AOROY Fe seTbhuy MOTA uoTFeTAed °7G oInbt gy 


o 


SSJIONY NOTIBIASZO 





62 02 "Gt ‘Ol °S “oO 
' . "9 
"O} 
a «) 
D 
cle 
es 
Dp 
~ 
oO 
wn 
+ 
. D 
nl z 
Cm 
Mm 
ra 
q> 
a 
“ato et 
wn 
"91 
CINWA OFIZIGIW 
GINA G3AH3SAD @ 
SS°IS) = MO19 SSHM 
N91S30 x04 = 033946 
‘02 


L3ILNG YOLOY 


eas 





OL ‘2ETUIT rOWeAS |YQ Qe AQTOOTOEA TeTxw “ES SANhHTA 


(eae 7a silt oe i) Os Ab Xb 


“609 OSs “acs “US? "Con “Ost “ore 

nS En NE eee '@ 
O01 

e 

© 

“et 

® 
m] 

e 
9! 
© 

“el 


SIN WA 031910300 V 
S3NIWA O3AN7IS80 © 
SS°1S1 © MOT4d SSUM 
wW3IS30 20 © 0334S 


‘02 


PEPIN stGlaE ol ILS 


3) 


3INHiSTO IWIiGbY 


([Sa-IN]) 





“001 





‘RZOTUT 709e89S 3e SOTbuy MOTY 


“08 


SINWA OF19103Nd ¥ 
S3INIWA O3AN3SRhO © 
SS°ISt = MOTs SSUM 
W91S30 ZOZ = 033dS 


SJTONY MOTI 31N10S8N 
°09 


1a) Tan 


“OR 


HOULGILS 


aynTosqyv 


"02 


"pG eanbty 
‘O 

"9 

‘Ol 
me} 
»D 

‘zz 0 
aD 
2 
oO 
Wn 
—4 

i Db 

ho> 
i oe) 
0 
Zz 
9 
a & 

91 Aa] 
uw 

“eT 

‘02 


140 





MoOrsOUeehte | Ae1Gd Odeo Ss eVOrsIey oinssetay Peto 8166 ounoe. 


OLY JYNSSIUd WIiOdl 


ort sQ'l 00°} $6°0 06°0 





‘@ 
‘Or 
D 
a) 
2) = 
D 
T— 
| 
Le 2) 
4 
: D 
ht > 
™ 
m 
ae 
a) 
ae 
"9! mv 
Ww 
“@! 
SaNWA OF19103Nd ¥ 
Sanwa O3Au3S80 OD 
SS°1tSt = MO14 SSUW 
N9IS90 XOZ = 033905 
“02 


LAIN! YOLUIS 


141 





ubtsedq OL *’FeETINO 70ReFS SY Fe AXTOOTOSA TeTXW °9G aainbt 7 


oO 


(23S/14) LLJIOIZA IWIXY 





“009 “OSS “00S “OSH “Ooh "OSE "OO¢ m 
‘ot 
a 
D 
21 2 
2D 
Uae 
oO 
w 
4 
nm 2 
oO 
m 
ae 
> 
mn 
91 os 
et 
SINWA O3131034d V 
S3NWA OFAHISA0 D 
SS*tSt = MOTs SStH 
W9IS9O ZOOL = 0394S 
‘p2 


IPS! SOS eS 


142 


ee a 


== 





ubtsaqd 


° 


°° 


0 


L '3eTIANO 104e4S ye sEThuy MOTT syNnTosqy 


“oe 


SSTONY HOTA 31NTOSeY 


“02 “Oo! “0 


“OI- 


"Cee 


"LG 9anbty 


“OR - 





SINIWA OIID103d F 
SINIWA 0OFAH3SE0 O 
SSIS) = MOTs SSUR 
N9IG30 XO © 0339d¢ 


TBE shel al LS 


“01 


“2 


“hl 


“ol 


“et 


"02 


JSNBLSIO (BI OeY 


SHINT) 


t 
4 


= 


243 


<= 





ubtseq 304 ‘3eETaANO AOReIS Je saTbuy MOTY UuoTIeTASqd °8G BINH y 


°o 


S3JTONY NOTIVIASO 


‘$2 "02 “St ‘oO! Ss 6% 


“O1 


2} 


NBISTIO Widsyy 


“hl 


9 
to en) 
ss 


a 
(SSHINI) 


‘ai 


SINWA O319103Hd 

S3INWA 03AH3S80 @ 

SS°1S1 = MOVs SSUR ‘ 
W9ISIO ZOL = 033946 


L3TLNOD YOIBIS 


"02 


144 





ubtsoaq 


be-<al 


"009 


0 


8 


‘qoOTUI AOROY 3y. Fe ARTOOTSA TeTXY 


‘OSS 


SINTWA 0313103Ud 
S3N WA G3AH9Se0 © 
HS'AZt = MOV SSHW 
W91S30 208 = 03346 


(23S/14) ALIIOIZA IWIKY 


*COS “OSf 


LI INI YOLOY 


“O0n 


"OSE 





7G Garo ulonea) 


“006 


"OT 


“et 


nhl 


"el 


“02 


TO TOBY 


BiSstd 


JON 


(SSHING) 


145 


- 





ubtseq 308 “‘FETUL AOROU 7e soTbuy MOTY SATIeTSEN °09 eANbtTY 


S3VINU MOS JATIU VAY 


"00! “oe "09 “On “02 0 





‘@ 
‘Ol 
D 
2D 
“et Ss 
2D 
[7 
S 
Ww 
—4 
“\l 
CY 
my 
> 
0? 
aie 
°O1 mm 
wn 
“91 
CINIWA OF1310Id FY 
SINIWA O3ANISEO @ 
AS°HRL1 = MOV Suu 
N9IS3O 208 = 0334S 
02 


PEPIN eI BBs 


146 








ubtsed 208 ‘10304 eYF JO OTReY SANSSerg Te}IOL “19 sANHTA 


Oo 


OlLBY SYUNSS3Hd WI0L 


| S*1 n't €°l 2° ae 





‘at 
DD 
D 
‘21 = 
D 
— 
oO 
uw 
=4 
Fi D 
fil = 
© 
m 
Pa 
i an) 
=e 
‘St mt 
v 
‘ot 
S3INTHA 0313)03Ud V 
SINTBA GIAHISAD O 
hS'hit = MOTs ESUR 
N9IS3IO 208 = 63306 
‘a2 


LEPINE eh BN IEDs 


147 





ubtseq 30g ‘10304 ey JO ADUSTOTFJqa OTAeqeTpW °7Z9 aainbt 4 


°o 


HOLOU “LINIIIISS3 JILLYeyNIOY 





SSNTWA O313103Ud ¥ 
Sanwa O3Au3Sa0 @ 
NS *h2! = MOT14 SSHH 
w91S39 208 © 033dS 





ibe NP eSB Bhs 


“él 


“hI 


‘ay 


*e1 


“02 


3SNeLSTO TWIidkY 


(SSHINT) 


148 





ubtsed $08 ‘3eTINO AOVOY 3yF 7k AQTOOTOA TeTXY 


"009 


“OSS 


(93S/14) 


“00S 


ALIJOVWIA WIXY 


“OSh 


“00h 


“ose 


69 Oana 


“006 





SIN WA O313103Hd ¥ 
S3INWA OFAHISEO D 
RS*HLT = MOTs SSUW 
NS3I630 X08 = 033dS 


Pea 


YOLOY 


“Ot 


| 


“ht 


“91 


“et 


ity 


JINYLSIO WWldey 


(S3HINT) 


Lao 





ubtsoqd 308 ‘3eTINO 42070H Fe soTbuy MOTI eaT3zeTSxY ° PY BInPT A 


o) 


SAJTONWN MOTTA JATLYISY 


‘001 “oe “09 “On nO2 0 





‘01 
D 
D 
‘21.0 
xD 
2 
S) 
he 
wn 
~4 
mos 
mM 
m 
=z 
om 
= 
‘a1 mM 
Wy 
Ol 
SINWA 03191030 ¥ 
GaN WA O2ANIGAD O 
NSRLE © HOTS SCUM 
916320 X08 © 0334S 
“02 


L33LNO YOLOY 


50 





ubtsoq 308 ‘39T3NO AOZOY We saTbhuy MOTTA saAnTosqw °s9 eANhTty 


SJTSNY MOTI 3LNTOSEY 


"aot “098 ‘09 ‘on ‘O02 “0 

“9 
‘ot 

= >) 

D 

‘2t Oo 

2D 

t~ 

oO 

wn 

nar 

A D 

ni 2 

S 

mn 

a 

om | 

oe 

97-89 
‘el 





SIMA OFLFIGIE VY 
SINIWA O3AN3S80 © 
AS*HLZT = MOT SSHU 
N9OISSO 208 = 0339d¢ 


ib JO) veaaals 


sg) 





ubtseq 308 ‘3ETINO AORVOY 3e SOaTHuy MOTA vuoTReTASsq *99 ZANhTYA 


SITONY NOTIYIA3ZO 





"$2 ‘02 St ‘ol Ss ‘0 
"9 
“01 
wD 
D 
‘21 4 
wD 
| a 
oO 
es 
in 
—~4 
“hl z 
= 
my 
= 
=z 
9m 
= @ 
‘ear om 
w 
1 
SINIWA O313)03Nd 
SINIWA OF3AH39S80 @ 
AS‘RZ1 = MOVs SSuM 
N91S30 200 = 03346 
"02 


LAdNO 0104 


foa2 





ubtsaq 





ek 

“hh 

“el 

‘al 
CINTA OIIIIAW VF 
SINTWA O3AH39SA0 © 
NS'HLZE © MOT SSUW 
W91S30 208 © 0334S 

‘02 


LS INI YOLUIS 


2909 ‘}aeTUT 70R7eRS SYA We AQRTOOTSA TeTxXv °*L9 aInbHty 
(235/14) ALIIOVIA WIXY 
‘009 “OSS “00S “OSh "OO “ase “CoC 
ct ee eae ae te ee a ee 9 
“Ol 


NolSIG Wide" 


cle 


3Ni) 


(SSH 


ubtTseqd 30g ‘FOSTULI WOReAS We SOTHuy MOTY sanTOSqu °*8g9 aanhbtg 


SJIONY MOVS 3LNIOSey 





‘001 06 “O9 “Oh *O2 ‘0 
“Q 
‘01 
m4) 
D 
‘21 0 
yD 
& 
oO 
w 
~4 
“fh = 
cy 
m 
<= 
O 
Ir 
‘Ol m™ 
ol 
*@1 
GINWA O3LI3I03Ud Y 
SIN WA OZAN3ISGO © 
HS'h2Z} = MOTs SSUW 
NOIGIA 208 = 0393d6 
“02 


IPB NESE ES 


154 





ubtsed 308 ‘FeTUL TOReAS Fe OTJEY SANsserd TeIOL °69 sinbtq 


Orley JYNSS3Hd Widl 





orl so'l on’ t S60 bell —oa 
“Ol 
a 
D 
meee Se 
D 
i- 
(cs 
Ww 
=A 
: D 
hoo 
O 
hal 
=< 
ty 
2c 
(oie 
wn 
91 
SINWA O3139)03Hd ¥ 
SIN WA O2AvasSaO @ 
AS‘RZT © MOTI SSUK 
w99S30 200 = 03346 
‘02 


LIINI YOIBIS 


LS) 


- a 





ubtsoqd 


o 


(94$713! Stioe7an set xu 


‘009 “OSS “VOS “OSh “OOR "OSE “006 





SINWA 03131030d ¥ 
S3NIWA O3AHISEO © 
AS'het = MOTs SSHH 
W91S30 X08 = 033dS 


SPAY ei Islas 


‘ot 


“et 


“nt 


“Ol 


‘02 


208 '2ETINO JOReAS |YQ Qe AQTOOTAA TetIxw “OL eANbTY 


Wide" 


JNBLSTQ 


~ 
“dj 
a 


([SSHINT) 


E56 


=2)>-— ae 





ubtsaqd 


2) 


=< 


oO 


08 ’3eT3NO TOReAS We SeTbuy AMOTA 


“O€ 


SJ TONY MOTY JINIOSAY 


‘02 “ol “0 


SINWA 0319103Ud V 
SINWA O3AN3S5A0 @ 
AS'HLT &= MOTs SSUM 
N91S30 208 = 0334S 


OPE A MSNA eS: 


OY) ies 


aynTosqy 


“O2- 


“Eh Sant 


“Ob - 


‘O1 


“2t 


ht 


“02 


JINGLSIO IwlOby 


(SSHINT) 


Ty, 





og ‘2°T3NO AORWeRS 3e SEeTHhuy MOTY UOoTReTADG 


SJ TINY NOTIYIASZO 


“Se “02 “St "Ol "Ss 





S3NIA 0315103Ud 
SINHA O3AN3S80 @ 
RS'HLT = MOTs CSU 
W9ISIO X08 = 033dS 


ISIE te Els 


‘or 


Te4 


“hl 


“SI 


“6! 


“02 


‘CL OG bad 


INBLISIO WWIideY 


3 


(SSHINT) 


158 





abeyootd 36 ‘ubtseq 30g ‘3eTUI AOAVOY ey Ye ARTOOTOA TeTXWY ‘EL eAnbty 


(J3IS/14) ALTIOWA WIX 


°009 “OSS “00S “OS “OOh “OSE ‘00¢ | 


| cS ce oe eee Sees nee Gabe enti cee Gomes eee ees eee eee eee eee eee eee eee eee ite 9 
9 

8 ‘Oo! 

D 
e 
‘21 © 
is 
D 
f— 
© 
wn 
e ~_ 
fe Dw 
Lan Ga 
2 
m 
aa 
Fad 
© m 
rn 
91 in 
5 
6 
‘QT 


20°6 = YOIINs 39N9018 
SINWA O3131030d ¥ 
SINIWA 03AH36e0 @ 
WS°h2ZT = MOTs SSUH 
N91S3I 208 = 0334S 


LA INI YOLOY 


oo 








ObeyoOOTA &6 


‘ubtTsed 308 ‘3#eETINO AOROY ayW ye ARTOOTSA T[eRTXY 


“OSS “009 


eT Tg cena fern ona ren arma re oemanetnteny ceemanad 


10'6 = WO!IHs 294UNIOTE 
SINWA O319103uUd F 
SIN WA OFAN3SED © 
NS RET © MODY SSHW 
ND1C30 208 = 03306 


2257 2S a as Ay WI Xe 


“OSS ‘00S 


LaqiNO yOLOY 


“OSt 


“OOF 


“OS6 





"pL eanbtd 


‘0! 


“et 


“hl 


"02 


[OQ IWIideYy 


INT) YSNHIS 


(S 34° 


160 





OHEYOOTA $6 


nUDES OG 1306 ° donb tO1cso sua Je Aqnooion per ay 


ese “009 





Z0°6 © WOLIUS 35NNI0 
SIN WA 03191030d 


S3N WA O3Au3Se0 © 
WS'HLT © MOTs Sum 
NOIS3IO ZO98 © O33dS 


(33S/13) 


"ess 


IP SPIN 


LLYIOIGA WIXv 


YOLULS 


"GL 


“ot 


acu 


un 


61 


“@! 


‘a2 


O1InbHT gy 


JINBISTO WICH 


(SSHONT) 


ile) Ah 





abeyoota 36 ‘ubtsedq 308g “3eTINO AORVeAS SY 4e KATOOTaA Tetxy ‘9L eaANnbTY 


(J3S/13) LLEIO ISA 1IXY 


“002 “oss ‘009 ‘OSS ‘00S “Osh “ooh 

| SLT LT FR ETE, TL APC Ee Pe oe [saan fo Psa cenah emi ea fee Gy pen) ‘6 

‘oO! 
bs D 
i=) D 
‘21 oO 
| 
D 
iv 
e << 
wn 
—4 
ni 
Mm 
9 my 
z 
m 
2G 
© "at mm" 
uw 

@ 
@ 
“OI 


10°6 © YOIIU4 39¥N9078 
SINWA 031391030d ¥ 
SINIWA 0FAN3S90 © 
NS'hZ1 = MOTs SSUW 
M9OIS3IO 208 = 03346 


ee nU MSP cieelerls 


‘02 


ee 





APPENDIX A 


STORAGE ALLOCATION FOR THE PROGRAM MESHGEN 


The following list of pointers and variable names indi- 


cates the storage location of the first value of the corres- 


@emaing array: 


REAL 8 Variables 


Poameer 
Name 


op! 


J2 


J3 


J4 


J5 
J6 
J7 
J8 


JY 


(Array R8) 

Array Array 

Name Contents 

XSE Z Coordinates of the Super 
Elements 

YSe R Coordinates of the Super 
Elements 

ANSS Temporary Storage of the Z 


Coordinates of the Super 
Element Subdivision 


Wise Temporary Storage of the R 
Coordinates of the Super 
Element Subdivision 


ZC Z Coordinates 
RC R Coordinates 
PSI Stream Function 
B Blockage Factor 


End of Stored Values 
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REAL 4 Variables (Array 04) 


Pointer Array Array 
Name Name Contents 
M1 OHZC Z Coordinates of Streamwise 


Element Boundaries 


M2 OHRC R Coordinates of Streamwise 
Element Boundaries 


M3 OYVZC Z Coordinates of Transverse 
Element Boundaries 


M4 OVRC R Coordinates of Transverse 
Element Boundaries 


M5 OHZC1 Z Coordinates of Streamwise 
Element Boundaries 


M6 OHRC1 R Coordinates of Streamwise 
Element Boundaries 


M7 OLE Real 4 Z Coordinates of Mesh 
M8 ORC Real 4 R Coordinates of Mesh 
M9 End of Stored Values 


INTEGER 4 Variables (Array I2) 


Pointer Array Array 
Name Name COnECtEsS 
Ll NSEC Number of Columns per Super 
Element 
L2 NISE super Element Type 
L3 NTE Element Type Identifier 
L4 NODE GConmmeceivi ty Matrix 
L5 NBC Boundary Condition Nodes 
Lo End of Stored Values 
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APPENDIX B 


Prot iING OF THE OUTPUL VARIABLES FOR THE PROGRAM MESHGEN 


1 3¥5) 


FORTRAN Name Variable File 

ZC Z Coordinate TMESH DATA 

Re R Coordinate TMESH DATA 

B Blockage Factor TMESH DATA 

sass Inlet Static FLUID DATA 
PEesomine 

PTOT timlee Total FLUID DATA 
Pressure 

TEMP Inlet Static FLUID DATA 
Temperature 

VOT Inlet Total ELOLD pars 
Temperature 

RHOSTA Pibet otatac PLUD pare 
Density 

RHOTT tolet feta 1 FLUID DATA 
Density 

WDOT Inlet Mass Flow FLUID DATA 

CP Specific Heat FLUID DATA 

R Gas Constant FLUID DATA 

G FLUID DATA 

RPM RPM FLUID DATA 

VZI Inlet Axial PhULD DATA 
Velocity 

NODE Connectivity Matrix NODE DATA 

NTE Element Type ID NODE DATA 





NBC Boundary Nodes BC DATA 
Where wy Specified 


eo 1 Stream Function STREAM DATA 
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APPENDIX C 


STORAGE ALLOCATION FOR THE PROGRAM TURBO 


The following list of pointers and variable names indi- 
cates the storage location of the first value of the corres- 


Bomaing array: 


REAL 8 Variables (Array R8) 


Pointer Array Array 

Name Name Contents 

NP1l ZG Z Coordinates 

NP2 RC R Coordinates 

NP3 B Blockage Factors 

NP4 ALP Absolute Flow Angles 

NP5 BE Relative Flow Angles 

NP6 H Total Enthalpy 

NP? Hs Stature Hite tadipy 

NP3 VZ Ax1lal Velocity 

NP9 VR Radial Velocity 

NP10 VU Absolute Tangential Velocity 
NPil WU Relative Tangential Velocity 
NP12 Pot Current Streamfunction Values 
NP13 PSiO Previous Streamfunction Value 
NP14 F Right-hand Side Vector 

NP15 Ris Memgonary Storage ;£or F 

NP16 RHO Static Density 
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NP1?7 
NP18 
NP19 
NP20 
NP21 
NIE 2 2 
NP23 
NP24 
NP25 
NP26 
NE27 
NP28 
NUBIA, 


NP30 


REAL 4 Variables 


Pointer 
Name 


NPO1 


NPO2 


NPO3 


NPO4 


NPOS 


RHON 


WRL 


ETA 


EM 


DEV 1 


PRAT 


TEMP 


TOT 


PRE oo 


aoa 


ROL 


HR 


ENTROP 


(Array 04) 


Array 
Name 


OVEL 
ORC 
OBE 


OALP 


Not Used 

Angular Momentum Vector 
Adiabatic Efficiencies 
Stiffness Matrix Elements 
Rotor/Stator Deviation Angles 
Total-to-Total Pressure Ratio 
Static Temperature 

Total Temperature 

Static Pressure 

Total Pressure 

Total Density 

Rhothalpy 

Entropy 


End of Stored Values 


Array 
Contents 


Velocity at a Given Station 
R Coordinates 

Relative Flow Angles 
Absolute Flow Angles 


End of Stored Values 
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INTEGER 4 Variables 


Pointer 


Name 


Nett 


NE ie 


NPI3 


NPI4 


NPIS 


Array 


Name 


NFS 


NBC 


NTE 


NODE 


(Ammay i2) 


Array 
Contents 


Nodes for F Calculation 
Boundary Nodes 

Element Type Identifier 
Connectivity. Matrix 


End of Stored Values 
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APPENDIX D 


meee iNG OF THE OUTPUT VARIABLES FOR THE PROGRAM TURBO 


Listing Name 


FL 


VZ 


VR 


R 


DENSITY 


WRL 


gill 


VOM 


WT 


HS 


ieMP 


cee). 


PRESS 


PLOT 


RHOT 


ALPHA 


BETA 


Variable 
Stream Function 
Axial Velocity 
Radial Velocity 
Radius 
Static Density 
Angular Momentum 
Total Enthalpy 


Absolute Tangential 
Veroe uEy 


Relative Tangential 
Velocity 


Static Enthaipy 
Static Temperature 
Total Temperature 
Static Pressure 
Total Pressure 
Total Density 
Absolute Flow Angle 


Relative Flow Angle 
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lbm/sec 
ft/sec 
ft/sec 
inches 
lbm/f£t> 
£t*/sec 
BTU/1bm 


Et/sec 
ft/sec 


BTU/1lbm 
mine 
rps 
psia 
psia 

3 
lpm, fe 
Degrees 


Degrees 
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te eee 


MESHCEN FCRTRAN AL NAVAL PESTGPAQLATE SCHOOL 
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FEEH EERE ASHE HEE EHH ASH EEE EE EEE EHR SEHR EERE HERE EH E HE HEHEHE HHH He tH Hee ett 
PEE EEE AEE EEE EH ELH EEE EEE HEH EE EEE HEHE HEE EE HEHEHE FEE EEE HEHEHE Ee Hee ee eee 4 


++ ++ 
++ PRCGRAM MESFCE t+ 
++ ++ 


HEATH HAE EAE EE EE EEA EEE EH HEE EHED EEE EEE HEHEHE HEH HE HEHE HEHE HH Eee etete ett 
HE EE EE EEE HE EEE EE THEE EEE EEE EE HEE HEHE EHHE FE HAASE ERE TE EEE EEE HH HE Her et 
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HN + + 
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++ ++ 


SHEVA DEAE EEA EA ERE THEE FE HEE EE EE HH AE EHH EEE EH EE EEE EH EEE EH HEEE HEHEHE AHHH Ht 
SHAPE EEE HE AEE EHR HEHEHE HEHEHE EE HEHEHE EE HEHEHE HEHEHE FEE EEE EEE HEHE EHH EH He te Hose 


IMPLICIT REAL FECA RWI, P HZ), RHAL #404) 
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fees MESFGEN FERTFER A] NAVAL POSTGRADLATE SCHOOL 


20C Pte eeks PENG RY SPHC= AVETLABLE = REAL =",15-2X,' INTEGER = * 
Lek 7-77) 
FETURN 
EAC 
C 
: 
See Ue cee Sl 2 oe y cc eX A SEs YYS So ZC yekC NSEC, RTSE NTE, 
PeaASE27 eA SEN2, MFI, FANS) 
C FEE AE AEEH EF EEE FE HAE HEHEHE ER EERE EE FEE HEH HEHEHE FEE EE EE HEHEHE Hee He eee tee t+ 
C FHP HEHE AEA HEHEHE HH EHEH EE EHEHHE HEHE HEHEHE HEHE EH EEE HEHE EEE HEHEHE H HEH HHH HHH HHH HH HH HH + 
C ++ ++ 
C ++ TEIS SULEFCUTIAES FRCTICES A RECTANGULAR EITCHT-NSDED ++ 
= ++ ISECPASP SMe T41C GPin CE ARPSITPAFY DIAFASIO' 4X N GIVEN ++ 
C ++ ie OC mieeo 9 Nth es Se Tree ESE’S "SUPER SLEWENTS* + + 
C +4 ++ 
C SHH + SAFE PH HE HEHE eH eH HH HEHEHE HEE HEHEHE HHH HHH HH FASE SEH HH HHH HHH HH He HH He eee eet 
C FHF HES EH EH HEATER FEA HEE HEE EE EEE EFHE HEHE HEHE HEHE HEHEHE EE HEE HEHEHE HEHEHE HE HHH H+ + 
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IMELICIT FEALFE(A-H,F-2) 
eGo oe ee th bei p lus CC» “RE Le licl » MN 1y YNZ, 4NG »YNG NSE2 gi2 
INTEGER *4 TANS 
GCOVVEtley iio 84 NR I GE eNC Ly IC 
COM4ON SMOCCUINT J ee On YR GI 
Mer went, PCE IN / oe bp COL, hey I2QTC»NSTAT 
es Oe ey oes Koei he YY SEC LI eZCIL) sRC(lL), 212) 
A Celie Ss (EC Nous EGU)» STSECL Py Mi e( 1) 
PeGl ano wecen) ClOIC 2C 
Ses i: 224744 C 
PeOOW2@ ie) 2C Cl) .FE (1) .a( 1) 
2 Pes se eS tt) 
rt CGhT ive 
Cette LOG 
z0 YC IV? = FLOAT(NRK 1) 
YeEeIV2 = FLEAT (#7) 
I = 1 
TJ = 1 
It = 1 
Iv = } 
ASTAT = C 
AFCTC = € 
SG tO Nw = t,hf 
Pelt) = UelC 
SO Ny) =O 6 tC 
10 CCATINGE 
COO ies eee eg a 
PVE S16 Os) Gua | 
Nee cel awe cot Jet) ) Loo yl 20s 140 
14C eter’ bee 1 Ge I) So 
ROTC. S11 
IM = € - i 
Gli: 126 
15 C Node. = 1k 
Itf4 = C€ - [M 
120 Iy = BSEC(S) 2 ‘ARC 
OG 1 ot Ss ie lt 
NTE( In) = *TSE(J) 
lle Le oa 
11¢€ CONTIPtLE 
Neal} =12 4 NSCC(J) + 1 
NG. = NCE] 
NE = NG) ] 
MEIVYL = FPEGCAT (PEC ]) 
XOTVE = FLOCATOSC) 
Mee Se) 
Kee NSE 1) 
Ke = PSqC1 4 2) 
MG = X¥SF( [1 4 2) 
Yl = YSe?t) 
YO .=s Ytrl i + 1) 
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Vise yottt 4° 3) 
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YCIF2 = yl - Y32 
YOIF2 = Y2 - Y4 
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>CC YCU bISrF TOC FNTFR TRE FLOT SEQUENCE? 
ee oY eee oe NC 

>OC YOU bITSKR TCO CCATINUE FLOT SEQUENCE? 
I= YES sc = NE. 

Pees YOU wISR FC COATINUE PLOT SCQUENCE? 
1 = YES; 2 = NC. 

Peeve alse TE CEATINUE FLOT SEQUENCE? 
Me = YES: <2 = KGe 

>O0C YOU wISH TC CLCATINUE FLOT SEQUENCE? 
=e YE Spo = Ls 

>CC YOU wISk TO CCATINUVE FLOT SEQUENCE? 
fees Ve Se = ACS 

PCC YCU KISH TC CCATINUE FLOT SEQUENCE? 
lL = Yc ¢ 2: NCe 

>CC YCU bISH TO CCATINUE PLOT SEQUENCE? 
1 = YES; 2 = NC. 

pa YS) VISE FE COAG TINUE PLOT SFOUVENCE? 
oss 2 = NC. 

DEG YCU wISk TC CCATIAUE PLOT SEQUENCE? 
Meee YES, ec = NCo 

pre YC) swalSh TC COAT INUE PLOT SEQUENCE? 
fe YESs 2 = NC, 

>COC YOY » ISH TG CCATINUE PLOT SFQUENC E? 

: 1 = YESs+ 2 = AC. 

meee yoU wilSH TC COATINJE PLOT SEQUENCE? 
Meee YES; . = NCo 

>CC YCU awISP TC CCATINUE PLOT SEQUENCE? 
Peo Yes, <c = AC. 


246 





exe DISLRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93940 


Department Chairman, Code 67 
Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 


Bececcor, Turbopropulsion Laboratory, 
Code 6/7Sf 

Naval Postgraduate School 

Monterey, California 93940 


LT J. Ferguson 

Naval Air Systems Command 
AIR-536 

Washington, D.C. 20361 


feed. Okiishi 

Department of Mechanical Engineering 
Iowa State University 

Ames, Iowa 50011 


247 


No. Copies 


20 











Thesis 





hlktees 3 


F25/23 Ferguson 


ooak 


Finite element pro- 
gram for calculating 
flows in turbomachines 
with results for NASA 
Task-1 compressor. 





